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■ Abstract 
Q-i. 

. We develop a new, mathematically precise framework for treating the effects of nonlinear phe- 

nomena occurring on small scales in general relativity. Our approach is an adaptation of Burnett's 
q ' formulation of the "shortwave approximation" , which we generalize to analyze the effects of mat- 

ter inhomogeneities as well as gravitational radiation. Our framework requires the metric to be 
close to a "background metric", but allows arbitrarily large stress-energy fluctuations on small 
scales. We prove that, within our framework, if the matter stress-energy tensor satisfies the weak 
energy condition (i.e., positivity of energy density in all frames), then the only effect that small 
. scale inhomogeneities can have on the dynamics of the background metric is to provide an "ef- 

fective stress-energy tensor" that is traceless and has positive energy density — corresponding to 
the presence of gravitational radiation. In particular, nonlinear effects produced by small scale 
inhomogeneities cannot mimic the effects of dark energy. We also develop "perturbation theory" 
off of the background metric. We derive an equation for the "long-wavelength part" of the lead- 
ing order deviation of the metric from the background metric, which contains the usual terms 
occurring in linearized perturbation theory plus additional contributions from the small-scale in- 
homogeneities. Under various assumptions concerning the absence of gravitational radiation and 
the non-relativistic behavior of the matter, we argue that the "short wavelength" deviations of the 
metric from the background metric near a point x should be accurately described by Newtonian 
gravity, taking into account only the matter lying within a "homogeneity lengthscale" of x. Finally, 
we argue that our framework should provide an accurate description of the actual universe. 
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I. INTRODUCTION 



It is generally believed that our universe is very well described on large scales by a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) model. However, on small scales, ex- 
tremely large departures of the mass density from FLRW models are commonly observed, 
e.g., on Eartffl, we have 5p/p ~ 10 30 . Nevertheless, common sense estimates |2j suggest 
that (a) the deviation of the metric (as opposed to mass density, which corresponds to 
second derivatives of the metric) from a FLRW metric is globally very small on all scales 
except in the immediate vicinity of strong field objects such as black holes and neutron 
stars, and (b) the terms in Einstein's equation that are nonlinear in the deviation of the 
metric from a FLRW metric are negligibly small as compared with the linear terms in the 
deviation from a FLRW metric except in the immediate vicinity of strong field objects. 
These common sense estimates together with the fact that the motion of matter relative 
to the rest frame of the cosmic microwave background is non-relativistic strongly suggest 
that (1) the large scale structure of the universe is well described by a FLRW metric, (2) 
when averaged on scales sufficiently large that \Sp/p\ <C 1 — i.e., scales of order 100 Mpc in 
the present universe — the deviations from a FLRW model are well described by ordinary 
FLRW linear perturbation theory, and (3) on smaller scales, the deviations from a FLRW 
model (or, for that matter, from Minkowski spacetime) are well described by Newtonian 
gravity — except, of course, in the immediate vicinity of strong field objects. 

The above assumptions underlie the standard cosmological model, which has been re- 
markably successful in accounting for essentially all cosmological phenomena. Thus, there is 
good empirical evidence that assumptions (l)-(3) are at least essentially correct. Neverthe- 
less, the situation is quite unsatisfactory from the perspective of having a mathematically 
consistent theory wherein the assumptions and approximations are justified in a systematic 
manner. Indeed, it is not even obvious that assumptions (l)-(3) are mathematically consis- 
tent js], J]. In particular, nonlinear effects play an essential role in Newtonian dynamics, e.g., 
the fact that the Earth orbits the Sun arises from Einstein's equation as a nonlinear effect 
in the deviation of the metric from flatness. It is clear that one would get an extremely poor 
description of small-scale structure in the universe if one neglected the nonlinear terms in 
Einstein's equation in the deviation of the metric from a FLRW model; for example, galaxies 



1 If we go to scales of atomic nuclei, then Sp/p ~ 10 
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would not be bound. But if one cannot neglect nonlinear terms in Einstein's equation on 
small scales, how can one justify neglecting them on large (i.e., ~ 100 Mpc or larger) scales? 
In addition, since it is not clear exactly what approximations are needed for assumptions 
(l)-(3) to be valid, it is far from clear as to how one could go about systematically improving 
these approximations. 

Indeed, it is far from obvious, a priori, that nonlinearities associated with small-scale 
inhomogeneities could not produce important effects on the large-scale dynamics of the 



FLRW model itself, as has been suggested by a number of authors \5-\i7\ as a possible 
way to account for the effects of "dark energy" without invoking a cosmological constant, 
a new source of matter, or a modification of Einstein's equation. In fact, the example of 
gravitational radiation of wavelength much less than the Hubble scale illustrates that it is 
possible, in principle, for small-scale inhomogeneities in the metric and curvature to affect 
large-scale dynamics. The dynamics of a FLRW model whose energy content is dominated 
by gravitational radiation will be very different from one with a similar matter content but no 
gravitational radiation. It is the nonlinear terms in Einstein's equation associated with the 
short-wavelength gravitational radiation that are responsible for producing this difference 
in the large-scale dynamics. Although common-sense estimates indicate that similar effects 
on large-scale dynamics should not be produced by nonlinear effects of small-scale matter 
inhomogeneities in our universe, it would be very useful to have a systematic and general 
approach that can determine exactly what effects small-scale inhomogeneities can and cannot 
produce on large-scale dynamics. 

The main approach that has been taken to investigate the effects of small-scale inhomo- 
geneities on large-scale dynamics has been to consider inhomogeneous models, take spatial 
averages to define corresponding FLRW quantities, and derive equations of motion for these 



FLRW quantities 
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19l | . Since, in particular, the spatial average of the square of a quan- 
tity does not equal the square of its spatial average, the effective FLRW dynamics of an 
inhomogeneous universe will differ from that of a homogeneous universe. However, a ma- 
jor difficulty with this approach is that, when the deviations of the metric from that of 
a FLRW background are not very small, it is not obvious how to interpret the averaged 
quantities in terms of observable quantities. For example, if the total volume of a spatial 
region is found to increase with time, this certainly does not imply that observers in this 
region will find that Hubble's law appears to be satisfied. Further serious difficulties with 



this approach arise from the fact that the notion of averaging is slicing dependent and the 
average of tensor quantities over a region in a non-flat spacetime is intrinsically ill defined. 
In addition, the equations for averaged quantities that have been derived to date are only 
a partial set of equations — they contain quantities whose evolution is not determined — so it 
is difficult to analyze what dynamical behavior of the averaged quantities is actually possi- 
ble. This difficulty is well illustrated by a recent paper of Buchert and Obadia |20|, where 
they suggest that inflationary dynamics may be possible in vacuum spacetimes. However, 
this conclusion is drawn by simply postulating that a particular functional relation holds 
between certain averaged quantities under dynamical evolution. In fact, Einstein's equation 
controls the dynamical evolution of these quantities — so one is not free to postulate addi- 
tional relations — but the restrictions imposed by Einstein's equation are not considered. 

The main purpose of this paper is to develop a framework that allows us to consider space- 
times where there can be significant inhomogeneity and nonlinear dynamics on small scales, 
yet the frameworlcl is capable of describing "average" large-scale behavior in a mathemati- 
cally precise manner. We seek a framework wherein the approximations are "controlled" in 
the sense that they can be shown to hold with arbitrarily good accuracy in some appropriate 
limit. The results obtained within this framework will thereby be theorems, and the only 
issue that can arise with regard to the applicability of these results to the physical universe 
is how close the physical universe is to the limiting behavior of the theorems, in which the 
results hold exactly. 

The situation that we wish to describe via our framework is one in which there is a 
"background spacetime metric", g^ b \ that is supposed to correspond to the metric "aver- 
aged" over small scale inhomogeneities. In the case of interest in cosmology, g^j would be 
taken to be a metric with FLRW symmetry, but our framework does not require this choice, 
and no restrictions will be placed upon g^ until section IV. The difference, h a f, = g a b — g^ , 
between the actual metric g a b and the background metric is assumed to be small everywhere. 
This precludes the existence of strong field objects such as black holes and neutron stars, 
but even if such objects are present, by replacing these objects with weak field objects of 
the same mass, our framework should give a good description of the universe except in the 



2l| (see also 



22]), but our 



2 Our framework will have some significant similarities to the approach of 
assumptions will be considerably more general and our results will have considerably wider applicability. 
Our assumptions also will be stated much more precisely. In addition, we will develop perturbation theory 
within our framework. 



immediate vicinity of these objects. However, even though our framework requires h ab to be 
small, derivatives of h ab (say, with respect to the derivative operator V a of the background 
metric g^J) are not assumed to be small. Specifically, quadratic products of V c h ab are al- 
lowed to be of the same order as the curvature of g^ b \ Thus, a priori, such terms are allowed 
to make a significant contribution to the dynamics of itself. Finally, no restrictions are 
placed upon second derivatives of h ab . In particular, if matter is present, the framework 
allows 5p/p^> 1. 

How can one formulate a mathematically precise framework where approximations such 
as the smallness of h ab = g a b—g$ are "controlled" in the sense that limits can be taken where 
they hold with arbitrarily good accuracy? The basic idea is to consider a one-parameter 
family of metrics g a &(A) that has appropriate limiting behavior as A — > 0. To illustrate this 
idea, consider the much simpler case of ordinary perturbation theory, wherein one wishes to 
describe a situation where not only is g a b — gfb small, but all of its spacetime derivatives are 
correspondingly small. To describe this in a precise way, we can consider a one-parameter 
family of metrics g a b(X,x) that is jointly smooth in the parameter A and the spacetime 
coordinates x. The limit as A — > of this family of metrics clearly exists and defines the 
background metric g^ b \x) = g a b(0,x). If we assume that g a b(X) satisfies Einstein's equation 
for all A > 0, it follows immediately that g^ b also satisfies Einstein's equation. The first 
order perturbation, 7^, of gj^ is defined to be the partial derivative of 5W>(A) with respect 
to A, evaluated at A = 0. It satisfies the linearized Einstein equation, which is derived by 
taking the partial derivative with respect to A of Einstein's equation for <? a &(A), evaluated 
at A = 0. More generally, the nth order perturbation, 7^\ of g^J is defined to be the nth 
partial derivative of g a &(A) with respect to A evaluated at A = 0, and the equation it satisfies 
is derived by taking the nth partial derivative with respect to A of Einstein's equation. The 
perturbative equations for the metric perturbation at each order hold rigorously and exactly. 
Of course, the issue remains as to how accurately an nth order Taylor series approximation 
in A describes a particular metric g a f>(A) for some small but finite value of A. This issue 
may not be easy to resolve in any specific case. Nevertheless, even if the accuracy of the 
Taylor approximation cannot be fully resolved, it is far more satisfactory mathematically to 
derive rigorous results for the perturbative quantities than to make crude arguments about 
g a b based on the assumption that h ab = g ab — g^ is "small." 

To obtain a mathematically precise framework that can be applied to describe situations 
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relevant for cosmology, we also wish to consider a one-parameter family of metrics g a b{X) 
that approaches a smooth background mt 
require that first spacetime derivatives of 



that approaches a smooth background metric g^} as A — > 0. However, we do not want to 



h ab (\) = g ab (\) - (1) 

go to zero as A — > 0. Indeed, in order to capture the effects we are interested in, it is 
essential that, a priori, the framework allow quadratic products of derivatives of h a b(X) to 
be of the same order as the curvature of g^ b in the limit as A — » 0. This suggests that 
we should consider a one-parameter family wherein, as A — > 0, the deviations of g a b(X) 
from g^J simultaneously become of smaller amplitude and shorter wavelength, in such a 
way that first spacetime derivatives of h a b(X) remain bounded but do not necessarily go to 
zero. If h ab (X) — > as A — > but spacetime derivatives of h ab {\) do not go to zero, then it 
is easy to see that spacetime derivatives of h a b(X) cannot converge pointwise (i.e., at fixed 
spacetime points) as A — > 0. However, spacetime derivatives of h a b(X) will automatically go 
to zero when suitably averaged over a spacetime region; more precisely, their "weak limit" 
exists and vanishes. Similarly, although we cannot require that quadratic products of first 
spacetime derivatives of h a b(X) approach a limit at fixed spacetime points as A — > 0, it is 
mathematically consistent to require that the weak limit of these quantities exists. As we 
shall see in the next section, a certain combination of weak limits of quadratic products of 
first spacetime derivatives of h a b(X) acts as an "effective stress energy tensor", which affects 
the dynamics of the background metric g^ h . In this way, the possible effects on FLRW 
dynamics of small-scale inhomogeneities — which are required to be of small amplitude in 
the metric but may be of unbounded amplitude in the mass density — can be studied in a 
mathematically precise manner. 

In fact, the issues we confront in attempting to treat the effects of small-scale mass 
density fluctuations in cosmology are very similar to the issues arising when one attempts 
to treat the self-gravitating effects of short- wavelength gravitational radiation. In the latter 
case, one is interested in considering a situation where the amplitude of the gravitational 
radiation relative to some background metric g^ b is small, but the "effective stress-energy 
tensor" of the gravitational radiation — i.e., products of first spacetime derivatives of (g a b — 
^ab) — ^ S comparable to the curvature of g ® b . A "shortwave approximation" formalism was 



developed by Isaacson 
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24| (see also pp. 964-966 of [25|) to treat this situation. The 



shortwave approximation was put on a rigorous mathematical footing by Burnett 26] , who 
derived the equations satisfied by g^} by considering a one-parameter family of metrics g a fc(A) 
with suitable limiting behavior. In this paper, we shall generalize Burnett's formulation of 
the shortwave approximation by allowing for the presence of a nonvanishing matter stress- 
energy tensojfl T a b. By following Burnett's approach, we shall derive an equation for the 
"background metric" , , which takes the form of Einstein's equation with an "averaged" 
matter stress-energy tensor and an additional "effective stress-energy" contribution arising 
from the small-scale inhomogeneities. 

One of the main results of our paper, proven in the analysis of section II, is that if the 
true matter stress-energy tensor T ao satisfies the weak energy condition (i.e., if the energy 
density is positive in all frames), then the effective stress-energy tensor appearing in the 
equation for g^ must be traceless and must have positive energy density — just as in the 
vacuum casej. In other words, no new effects on large-scale dynamics can arise from small- 
scale matter inhomogeneities; the only effects that small-scale inhomogeneities of any kind 
can have on large-scale dynamics corresponds to having gravitational radiation present. Our 
analysis makes no assumptions of symmetries of the background metric g a °J and makes no 
assumptions about the matter stress-energy tensor T ao other than that it satisfies the weak 
energy condition. However, if g^ is assumed to have FLRW symmetry, then our results 
establish that, within our framework, the only effect that small-scale inhomogeneities can 
have on FLRW dynamics corresponds to the additional presence of an effective P = |p fluid 
with p > 0. In particular, within our framework, small-scale inhomogeneities cannot provide 
an effective source of "dark energy" . 

In cosmology, in addition to analyzing the dynamics of the background FLRW spacetime, 
one is interested in analyzing the deviations from the FLRW background. In section III, 
we undertake a general analysis of perturbation theory within our framework. It is not 
straightforward to do this because g a b{^) is not differentiable in A at A = 0, so there is no 
notion of a "linearized metric perturbation" in our framework. However, the weak limit as 
A —7- of [g a b{ty — 9 a °b]/^ ma y ex ist, and, under the assumption that it does, this limit 
defines a quantity which corresponds closely to what is called the "long wavelength 
A generalization of the Isaacson formalism to include matter under Newtonian assumptions was previously 
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considered by Noonan 

4 Our result on the positivity of energy densit y w ithin this framework is new, i.e., it was not previously 



shown to hold in the vacuum case by Burnett 
case was shown by Isaacson 



26] . Positivity of the effective energy density in the vacuum 



24] only under an ad&itional WKB ansatz. 



part" of the metric perturbation in other analyses (see, e.g., [28]). We also write 

^i?(A) = MA)-A 7 2 ) (2) 

( S) 

and refer to h ab as the "short wavelength part" of the deviation of the metric from a FLRW 
model. (Note that in our framework, h^ b and 7^5 have precise mathematical definitions.) 
Our goal is to derive the equations satisfied by 7^ as well as to determine (A) to accuracy 
0(A). This will yield the spacetime metric to accuracy 0(A). 

In section III, we will systematically derive the equations satisfied by 7^ and ft^f (A) in 
a completely general context. By taking the weak limit of 1/A times the difference between 
the exact Einstein equation for g a b(X) and the effective Einstein equation for g^ b , we ob- 
tain an equation for 7^ corresponding to that arising in ordinary linearized perturbation 

theory. However, in addition to the familiar terms appearing in ordinary linearized per- 

( S) 

turbation theory, this equation contains additional "source terms" arising from h ab and, if 
gravitational radiation is present in the background spacetime, this equation also contains 
additional terms linear in 7^ and quadratic in We also obtain additional relations 

between quantities appearing in this equation by taking weak limits of Einstein's equation 
multiplied by h ab /A and by h ab h cd /A. These relations are used to simplify the perturbation 
equation for 7^ . Finally, we write down Einstein's equation for h^ b (A). In the vacuum 
case, we consider the simplifications that can be made to this equation if one is interested 
only in determining h^ j) (X) to sufficient accuracy to obtain qffi . We compare our approach 



to that of Isaacson 
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24j and subsequent works (see, e.g., [251]). 
In section IV, we apply our general perturbative analysis to cosmology. We introduce the 
"generalized wave map gauge" in subsection IVA. In subsection IVB, we make additional 
assumptions concerning initial conditions and the Newtonian nature of the deviation, 5T a t = 
T a b{X) — T^ b \ of the stress-energy tensor from a FLRW model. We argue that, for small 
A, h ab (A) should be well approximated (in the wave map gauge) by the Newtonian gravity 
solution, whereby one needs only take into account the "nearby" matter. In contrast to 
the rest of the paper — where all of the assumptions are stated in a mathematically precise 
manner, and all of the results are theorems — in subsection IVB we provide only a sketch 
of the assumptions needed, and many of our arguments have the character of plausibility 
arguments rather than proofs. In subsection IVC, we simplify the equation for 7^ derived 
in section III by using our Newtonian assumptions. We show that this equation reduces to 
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the ordinary cosmological perturbation equation with an addi 
from h ab (A), in agreement with a result recently obtained by 



ional effective source arisinc 



In summary, in this paper we introduce a new framework for treating spacetimes whose 
metric is close to that of a background metric gj^ but is such that nonlinear departures 
from are dynamically important on small scales. We proceed by introducing a suitable 
one-parameter family of metrics g a f>(A) and deriving results in the limit as A — > 0. We prove 
that the small-scale inhomogeneities cannot affect the dynamics of the background metric 
gjfi except by the addition of an effective stress-energy tensor with positive energy density 
and vanishing trace, which can be interpreted as arising from gravitational radiation. We 
derive an equation for the "long-wavelength part" of the leading order deviation of the metric 
from g^ h , which contains the usual terms occurring in linearized perturbation theory plus 
additional contributions from the small-scale inhomogeneities. Finally, we argue that the 
small-scale deviations of the metric from gjj£ should be accurately described by Newtonian 
gravity. 

Of course, the real universe is not the limit as A becomes arbitrarily small of the type 
of a one-parameter family of metrics g a ft(A) considered here. Thus, our results do not ap- 
ply exactly to the real universe — any more than the results of an analysis using ordinary 
linearized perturbation theory would apply to a real situation. However, in section V we 
will argue that since the scales at which nonlinear dynamics are important in the present 
universe (i.e., scales much less than ~ 100 Mpc) are much smaller than the scale of the back- 
ground curvature (i.e., the Hubble radius ~ 3 Gpc), it seems reasonable to expect that the 
real universe will be accurately described by a "small A" approximation to <? Q b(A) within our 
formalism. We believe that our analysis thereby goes a long way toward providing a mathe- 
matically sound framework that can be used to justify the assumptions and approximations 
used in cosmology. At the end of section V, we will discuss how these approximations can 
be improved. 



Our notation and sign conventions follow that of [29] . Lower case Latin indices from early 
in the alphabet (a, b, c, . . . ) denote abstract spacetime indices. Greek indices denote com- 
ponents of tensors. Latin indices from mid-alphabet k, . . . ) denote spatial components 
of tensors. 
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II. DYNAMICS OF THE BACKGROUND METRIC 



In this section, we will give a precise statement of the assumptions that underlie our 
framework. We will then analyze the dynamics of the background metric g^J and prove that 
if the matter stress-energy tensor, T a i, satisfies the weak energy condition, then the "effective 
stress-energy" contributed by small-scale inhomogeneities must have positive energy density 
and vanishing trace. 

As explained in the Introduction, we wish to consider a situation wherein we have a one- 
parameter family of metrics g a bW that approaches a background metric g^J , but spacetime 
derivatives of g a fc(A) do not approach the corresponding spacetime derivatives of g^J . An 
example of the type of behavior that we have in mind is for components of h a b(X) = g a b{^) — 
g^ to behave like Asin(a;/A). In this situation, if we let V a denote the derivative operator 
associated with g^ b , we cannot have V c h ab (X) — > pointwise as A — > 0. However, suitable 
spacetime averages of V c /i a &(A) will go to zero. More precisely, if f cab is any smooth tensor 
field of compact support, we have 



provided only that h a h(X) — > locally in L 1 , where the volume element in this integral is 
that associated with g^ b \ If ([3]) holds for all "test" (i.e., smooth and compact support) 
tensor fields, f cab , we say that V c h a b(X) —> weakly. More generally, if A ai ... an (A) is a one- 
parameter family of tensor fields defined for A > 0, we say that A ai „ an (X) converges weakly 
to -Bai...a n as A — » if for all smooth / ai - a » of compact support, we have 



Roughly speaking, the weak limit performs a local spacetime average of A ai „ Mn (\) before 
letting A -> 0. 

As noted above, if g a b{X) converges to g^ b in a suitably strong sense — locally in L l 
suffices — then spacetime derivatives of h a b(X) automatically converge weakly to zero. How- 
ever, there is no reason why products of spacetime derivatives of h a b(\) must converge weakly 
at all, and, if they do converge, one would not expect them to converge to zero. (This latter 
observation is closely related to the fact that averages of products of quantities are not nor- 
mally equal to the product of their averages.) As discussed in the Introduction, we wish to 




as A -> 



(3) 
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consider a situation where first spacetime derivatives of h ab (X) remain bounded but do not 
necessarily go to zero as A — > 0. However, in order to have well defined averaged behavior 
in the limit as A — >■ we want the weak limit of the nonlinear terms in Einstein's equation 
to exist. As we shall see below, this will be the case if the weak limit of quadratic products 
of first spacetime derivatives of h ab (X) exists. 

All of the above considerations lead us to consider a one-parameter family of metrics 
<7a&(A), defined for all A > 0, satisfying the following conditions, which are straightforward 
generalizations to the non-vacuum case of the conditions considered by Burnett 26[ in his 
formulation of the shortwave approximation. In these conditions, V a denotes an arbitrary 
fixed (i.e., A-independent) derivative operator on the spacetime manifold M. For convenience 
in stating these conditions, we choose an arbitrary Riemannian metric e ab on M and for any 
tensor field £ Ql ... an on M we define |t ai ... a J 2 = e aibl . . . e anb H ai ,„ a J bl ... K . 

(i) For all A > 0, we have 

G ab (g(X)) + Ag ab (X) = 8nT ab (X) , (5) 

where T ab (X) satisfies the weak energy condition, i.e., for all A > we have 
T a b{X)t a (X)t b (X) > for all vectors t a (X) that are timelike with respect to g a b{X). 

(ii) There exists a smooth positive function C\{x) on M such that 

\h ah (\x)\ < Xd(x) , (6) 
where h ab (X, x) = g a b(X, x) - g ab (0, x). 

(iii) There exists a smooth positive function C*2(x) on M such that 

|V c MA,x)| <C 2 (x) . (7) 

(iv) There exists a smooth tensor field [i a bcdef on M such that 

w-lim [V a h cd (X)V b h ef (X)} = fi a bcdef , (8) 

A — ^0 

where "w-lim" denotes the weak limit. 

It follows immediately that fi a b(cd)( e f) = fiabcdef and fi a bcdef = Vbaefcd, and it is not difficult 
to show |26| that fi( a b)cdef = fiabcdef- It also is not difficult to see that if g a b(X) satisfies the 
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above conditions for any choice of V a and e ab} then it satisfies these conditions for all choices 
of V a and e ab . In our calculations, it will be convenient to choose V a to be the derivative 
operator associated with the background metric g a ° b = g a b(0), and in the following, we shall 
make this choice. We shall also raise and lower indices with gj^. 

As discussed in the Introduction, the key idea is that our actual spacetime, with all of its 
inhomogeneities, is described by an element of such a one-parameter family, at some small 
but finite value of A. By analyzing the limiting behavior of such one-parameter families at 
small A, we hope to attain an excellent approximate description of our universe. However, 
unlike ordinary perturbative analyses, our one-parameter family g a b{X) is not differentiable 
in A at A = 0, so we cannot define perturbative quantities or obtain useful equations by 
differentiation with respect to A. 

Our first task is to derive an equation satisfied by the background metric g a °J = g a b(0). 
This equation will follow directly from Einstein's equation (J5J) for g ab {X), using the general 
relationship between the Ricci curvature of g^j} and g ab (X), namely 

Rab(g {0) ) = Rab(g(X)) + 2V [a C c c]b - 2C d b[a C c c]d , (9) 

where 

C c ab = ^g cd (X) {V«0 M (A) + V b9ad {X) - V d g ab {\)} (10) 

and, again, we remind the reader that V a denotes the derivative operator associated with 
g a ° b \ so that V c g a ° b = 0. It follows immediately from that 

Rab(g {0) ) - lg ab (X)g cd (X)R cd (g i0) ) + Agab(X) 

= G ab (g(X)) + Ag ab (X) + 2V [a C\ ]h - 2C f b[a C e e]f 

-g ab (X)g cd (X)V [c C e e]d + g ab (X)g cd (X)C f d[c C e e]f , (11) 

and invoking the Einstein equation for A > 0, 

Rab(g {0) ) - l -g ab {X)g cd {X)R cd {g^) + A^ ab (A) 
= 87rT afc (A) + 2V [a C\ ]b - 2C f b[a C° e]f 

-g ab (X)g cd (X)V [c C e e]d + g ab (X)g cd (X)C f d[c C e e]f . (12) 

We now take the weak limit of both sides of ( TT2|) as A — » 0. It is easy to see that the 
weak limit of the left side exists and is equal to G ab (g^) + Ag a ° b \ Aside from the term 
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8irT a b(\), the terms on the right side of ffl2|) all contain a total of precisely two derivatives 
of h ab (X). These terms can be classified into the following types: (a) terms linear in h a b(X), 
corresponding to the linearized Einstein operator acting on h a t>(\); (b) terms quadratic 
in h a b(X), corresponding to the second order Einstein operator acting on h a t>(\); and (c) 
terms cubic and higher order in h a t>(\). The weak limit of terms of type (a) vanish by the 
type of argument leading to ([3]). The terms of type (b) depend upon h a b(X) either in the 
form V a h cd (X)V bh e f(X) — which has weak limit fi a bcdef — or in the form /icd(A)V a V&/ie/(A). 
However, since 

h cd (\)V a V b h ef (\) = V a [h cd (X)V b h ef (X)} - V a h cd (X)V b h ef (X) (13) 

and the weak limit of V a [h cd (X) Vfc/i e /(A)] vanishes, we see that the weak limit of h cd (X) V a V&/ie/ 
also exists and is equal to —fi a bcdef- Finally, it is easily seen that the weak limit of all terms 
of type (c) vanish. 

Since the weak limit of all terms in f[T2l apart from T ab (X) exist, it follows that the 
weak limit of T a b(X) itself also must exist, without the necessity to impose any additional 
assumptions on our one-parameter family. We write 

Tl 0) ^w-limT ab (A), (14) 

A— »0 

and we may interpret as representing the matter stress-energy tensor averaged over 
small scale inhomogeneities. Since T a b(X) satisfies the weak energy condition for all A > 
and since g a b(X) converges uniformly (on compact sets) to g a ° b \ it is not difficult to show 
that also satisfies the weak energy condition, i.e., T^t a t b > for all timelike vectors t a 
with respect to g a ° b \ The weak limit of (1121) then takes the form 

Gabig^ + Agf^SnT^ + Snt^, (15) 

where the "effective gravitational stress-energy tensor" t a ° b arises from the weak limit of terms 
of type (b) above and can be expressed entirely in terms of fiabcdef- A lengthy calculation 
(see 



J) yields 

i -r i 

o_j-(°) _ /_,, c de c d e , n cd e \ , cd ,, c d 

OHL ab ~ g \ r c de c d e ' r 1 c de J ' acbd c^r" ca bd 

4--u cd - -n c d +-;y c d --n cd (16) 
i A^ab cd 2^ ( ab ) c d 4^ cab d 2 abcd ' 



Note that t^} corresponds to the "Isaacson av erag e" of the second order Einstein tensor of 



h a b(X). It can be shown to be gauge invariant [26] 
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Following Burnett |26j, we decompose Aabcdef into two tensors a abcd ef = IJ>[c\[ab]\d\ef an d 



flabcdef = H(abcd)ef, SO that 

Aabcdef = ~ ^(o! c (ab)def + Oi e ( ab )f cd ~ Oi e ( cd )f ab ) 

+ Aabcdef + Pabe fed ~ Pcde fab ■ (17) 

Then, we have 



87rta6 = a a b cd + 2" cdab-^ia \cd\b) ~ ~^ab |« cd e ~ 2 « cdej- C 18 J 



Note that the right hand side is independent of (3 abcde f. 

• i Kluliiiirv 4"tir*~\ l/oir nrnnorfiao r\T ■ 

(0) 



The remainder of this section will be devoted to establishing two key properties of 
namely, t^ a a = and t a ° b t a t b > for any timelike vector t a of the background metric g a °). 



To prove these results, we need to obtain further information about a a b c def from Einstein's 
equation ffT2l . To do so, it is convenient to re- write this equation in "Ricci form" as 

R ab (g®) - Ag ab (X) = 8nT ab (X) - 4irg ab (X)g cd (X)T cd (X) + 2V [a C c c]b - 2C d b[a C c c]d . (19) 

We now multiply this equation by h e f(X) and take the weak limit as A — > 0. The left side 
clearly goes to zero. The weak limit of h e f(X)C d b ^ a C c c ^ is also easily seen to vanish. On the 
other hand, we have 

w-lim/i e/ (A)V [a C c c]fe = -w-limV [a /i| e/ |(A)C c c]6 

A — A — 

= -^[b\[ac]\d]ef9 y 

= -Uacfef ■ (20) 

Thus, we obtain 



beef 



47rw-lim/ ie/ (A) (t^X) - lg ab (\)g cd (\)T cd (\)) . (21) 



In particular, the weak limit of the right side of this equation must exist. By similar 
arguments starting with Einstein's equation in the form (TP2]) . it also follows that the weak 
limit as A of h e f(X)T ab (X) exists. We write 

Kefab = w-hm h ef (X)T ab (X) . (22) 

We will now show that the right side of ( 12 ip vanishes. We first prove the following lemma. 
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Lemma. Let A{X) be a one-parameter family of smooth tensor fields (with indices sup- 
pressed) converging uniformly on compact sets to A(0), and let B(X) be a one-parameter 
family of non-negative smooth functions converging weakly to B(0). Then A(X)B(X) — > 
A(0)B(0) weakly as A 0. 

Proof. Let F be a test tensor field — i.e., a smooth tensor field of compact support — with 
index structure dual to that of A. Let / be a smooth, non-negative function of compact 
support with / = 1 on the support of F. Then F = fF. We have 

KmJ(A(X)B(X)-A(0)B(0))F 

= hmj [(A(X)-A(0))B(X) + A(0)(B(X)-B(0))]F, (23) 

where contraction of the indices of A and F is understood. The second term is zero because 
B(X) — > B(0) weakly and A(0)F is a test function. On the other hand, we have 

(A(X) - A(0))B(X)F 



< J \(A(X)-A(0))F\\B(X)f\ 

< sup \(A(X) - A(0))F\ [ B(X)f, (24) 

where we have used the facts that B(X) > and / > in the second line. Now let e > 0. 
Since A(X) — >■ A(0) uniformly as A — > on compact sets, there exists Ai > such that 
sup^gsupp^ |(A(A) — A(0))F\ < e for A < Ai. Similarly, since B(X) — > B(0) weakly as 
A -)■ 0, there exists A 2 > such that J(B(X) - B(Q))f < e for all A < A 2 . Thus, for all 
A < min(Ai, A 2 ), we have 

J (A(X) - A(0))B(X)F <e(j B{0)f + , (25) 

Thus the first term in (123 D must be zero as well. □ 



The vanishing of n e f a b (see (1221) ) is a direct consequence of this lemma. To see this, let 
t a be an arbitrary timelike vector field in the metric . Then, we have 

K efab t a t h = w-lim h e f (X)(T ab (X)t a t b ) (26) 

A- — ^0 

Since by condition (ii) on our one-parameter families, \g a b(X) — g^\ < XCi(x), on any fixed 
compact region, we can find a Ao such that t a is timelike in the metric g a b{X) for all A < Ao- 
Since T ab (X) satisfies the weak energy condition, the function B(X) = T ab (X)t a t b is non- 
negative for all A < Ao- We previously showed that T ab (X) (and hence B) converges weakly. 
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Assumption (ii) also directly tells us that A a b(X) = h a b(\) converges uniformly on compact 
sets. Thus, from the lemma, we immediately conclude that 



n efab t a t b = (27) 

for all timelike t a in the metric g^j . However, since K e f a b = Kef{ab), we have K e f a bt a t b = for 
all timelike t a if and only if n e f ab = 0, as we desired to show. 

A similar argument establishes that the second term on the right side of (121j) also vanishes. 
Thus, we obtain 

<bcef = • (28) 



Consequently, ( 1T8|) simplifies to 



8^1? = a a \ d cd . (29) 



Taking the trace of this equation and again using (1281) we obtain our first main result of this 
section: 

Theorem 1. Given a one-parameter family g a b{\) satisfying assumptions (i)-(iv) above, the 
effective stress energy tensor t^ appearing in equation (|T5|) for the background metric 
is traceless, 

t {0)a a = . (30) 

We now show that t ® b satisfies the weak energy condition. Let t a be a timelike unit vector 
field with respect to g^J . We wish to show that tf b h a t b > 0. It is convenient to choose an 
orthonormal basis of with t a as the timelike vector. We will use Greek letters /i 7 u,p,... 
to denote spacetime components in this basis and Latin letters i, j, k, . . . from mid-alphabet 
to denote spatial components. Then 

8vr/°W 6 — n ^ 
0/JL ab b 1 — "OAtO^ 

jk 

~ a 0j0k 



i jk 



1 



l\^ijk +VjM J -^Hk \, ( 31 ) 



where in the second line we used the antisymmetry of a a b c def in the first two and the second 

*(o) 

L ab 
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two indices, and in the third line we used fl28|) . Thus, we have expressed t^)t a t h entirely 



in terms of spatial components of Habcdef, which will be useful for taking advantage of the 
positive defmiteness of the spatial metric. 

Aside from the tensor symmetries that arise directly from its definition, the only restric- 
tions on fiabcdef that we have at our disposal come from ( 12 8p . There is only one equation 
that can be derived from (1281) that involves only spatial components of ^ a bcdef, namely^ 

= - a ,% • (32) 

This yields 

Vij \i = jki ■ (33) 
Using this relation, we may re-write (l3"Tj) as 

STrtS^ = \ - -V,,," 7 " + - n;/ k k } . (34) 

For the remainder of our argument, we will work in a small neighborhood of an arbitrary 
point P G M. We will work in Riemannian normal coordinates x about P adapted to our 
orthonormal basis. Let /J> be a one-parameter family of smooth, non-negative functions 
with support contained in a 5-ball centered at P such that 

J [f 5 P (x)] 2 V=gd 4 x = l. (35) 



An explicit choice of f P is 



f P (x) = ——F{x/S) , (36) 



where F is any smooth, non-negative function of compact support contained in a ball of 
radius 1 satisfying J F 2 d A x = 1, but there is no need to make this particular choice. Instead 
of working with h a b(X), we introduce the quantity 

^ ab (5,X) = f 5 P h ab (X). (37) 

Note that for any fixed 5 > and A > 0, ipab is smooth and of compact support, so, in 
particular, t\} ab and all of its spacetime derivatives are in L 2 . Furthermore, it follows directly 
from the properties of h ab [\) and /J> that all components of ip ab converge uniformly to as 



5 Equation (l28l) states that the weak limit of h e f times the linearized Ricci tensor vanishes, i.e., it has the 
character of the linearized vacuum Einstein equation off of flat spacetme. The linearized Hamiltonian 
constraint — i.e., the vanishing of the time-time component of the linearized Einstein tensor — is the only 
component of Einstein's equation that can be expressed entirely in terms of spatial derivatives of spatial 
components of the perturbed metric. Equation (|32|) corresponds to this equation. 
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A — > at fixed 5. Similarly, components of V c ijj a b are uniformly bounded in A and x as A — > 
at fixed 5. Since ip a b is of fixed compact support, it follows immediately that ||^ a b|U 2 — > 
and || V c ip a b || l 2 is uniformly bounded as A — > 0, where H^a&IU 2 = / |^a&| 2 - 

Since fi a bcdef is smooth, it is obvious from (1351) and the support properties and positivity 
of if 5 P ) 2 that 

VnuafrpiP) = lim / ^ a/37P (x)(/p) 2 v/^d 4 x. (38) 
<5-»o J 

On the other hand, since (/J.) 2 is a test function, from the definition (jSJ) of fi a bcdef, a t each 
fixed 5 we have 

y A t ^a/37p(a ; )(/p) 2 v / -^ ( i 4 x = lim / d p h Q p(\)d u h^ p ( y \)(f P ) 2 y/^gd 4 x 

= lim J d^ipapd^py/^g d 4 x . (39) 

Here, in the first line, we replaced the derivative operator V a associated with with 
the coordinate derivative operator d a associated with Riemannian normal coordinates at P, 
making use of the fact that the definition of [i a bcdef is independent of derivative operator. In 
the second line, we used d^ipa^ = fpd/j,h a p + h^d^fp and the fact that the resulting terms 
in ( 1391) with no derivatives on h a p vanish in the limit as A — > 0. Taking the limit of ( 1391) as 
5 —7- 0, we obtain 

Htiuafa P {P) = hmlim / d^apdvipyp d 4 x , (40) 

o— >0 X— >0 J 

where we have used the fact that \J—g = 1 at P. Note that it is critical in this equation 
that the limits be taken in the order specified. 
The corresponding formula for is 

too(P) = -r lim lim / ^x {d t ^ jk d^ k ~ 2d 3 4 k %^ k + 2d 3 ^d k ^ k - d^/d^} . 

62% S— >0 A— ^0 J 

(41) 

where, in this equation, indices are raised and lowered with the flat Euclidean spatial metric 
rjij = diag(l, 1, 1) corresponding to the spatial components of at P. The major advantage 
of (1411) is that we can apply usual Fourier transform techniques to evaluate the integral 
appearing on the right side of this equation prior to taking the limit. If we had not "localized" 
^afc(A) by multiplying it by f P , the Fourier transform of h ab (\) could have been ill defined 
and, even if it were well defined, it would contain global information about h a b(X) rather 
than local information about the behavior of h a b(X) near P. 
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Our strategy now will be to prove that the quantity 

lim f d A x {d % 4> jk &i/S k - 28^%^ + 28^8^-8^/8^} (42) 

A— >0 J 

can be expressed as a sum of terms which are either positive, or which converge to zero as 
5—7-0. Positivity of %/ then follows immediately. We proceed by taking Fourier transforms 
with respect to the spatial coordinates x only, with the convention 

4'(*> *0 = J^yjl j d'x^t, x)e- lk x . (43) 

As previously noted, ipj k and all of its derivatives are obviously in L 2 , since ipjk is smooth 
and of compact support. Since the Fourier transform is inner product preserving in L 2 , we 
have 

*oo OP) = lim lim / dt d3fc \hk% k W k ~ 2k i k$ k i ~f* + 2kjk k ^f^ - htfifL*) . 

(44) 

We may decompose ipij into its scalar, vector, and tensor parts as 

$ij(t, k) = cr(t, k)kikj - 2<pq i:j + 2k(iZj)(t, k) + %(t, k) , (45) 

where k l Zi = = k l Sij, and s\ = 0. Here is the projection orthogonal to k l of the 
Euclidean metric on Fourier transform space. Since the various terms on the right side of (145]) 
are orthogonal at each k, it follows immediately that, for example, |</i>(fc)| 2 < |^> 1J '(fe)T/>y(ft). 
Since ipij and all powers of k l times ipij are in L 2 , it follows immediately that tp and all 
powers of k % times (p are in L 2 . Thus, we can freely take Fourier transforms of tp and all 
powers of k % times (p. Furthermore, since the L 2 norm of ipij — and, hence, the L 2 norm of 
ipij — goes to zero as A — > 0, it follows immediately that the L 2 norm of <p — and hence the L 2 
norm of tp — also goes to zero as A — > 0. Similarly, the L 2 norm of dip must remain uniformly 
bounded as A — > 0. Similar results hold for the other terms appearing on the right side of 
(S5J). 

Substituting the decomposition ( 145]) in (jUJ) and using the fact that ipy is real (which 
implies that ip(t, k) = tp(t, — k)), we obtain 

t$(P) = — lim lim / dt d 3 fc {kitfsjk&K - Bkik'tpip} . (46) 
32n <5->o A->o J i ) 

Thus, we see that the "tensor part", s^, of ip^ makes a positive contribution to the effective 
gravitational energy density t^. This may be interpreted as saying that, at leading order 
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within this framework, gravitational radiation carries positive energy density. The scalar 
a and the vector part Zi do not contribute at all, as might be expected from the fact that 
these quantities should correspond to "pure gauge" . Finally, the scalar (p makes the negative 
contribution 

2E V = -^- J &t4?kk i k i <p$ (47) 

to the effective energy density. 

In order to interpret the meaning of (p and E v , we note that by f )45l) . <p satisfies 

Ak%<p = -tfhftj + titifaj . (48) 

In position space, (148]) becomes 

4d%ip = -d i d^ j j + d^ipij . (49) 

To put the right side of this equation in a more recognizable form, we return to Einstein's 
equation (|T2|) and consider its normal-normal component relative to a t — const surface. 
This corresponds to the Hamiltonian constraint equation, which has the property that the 
only terms containing second spacetime derivatives of h a b involve only spatial derivatives 
of spatial components of h a b- To obtain this equation from ( 1T21) . we raise both indices 
with g ab (X) and then take the 00 component. Since we are working in Riemannian normal 
coordinates about P we also express the background metric as 

gfl = Va? - \R a ^x a x p + 0(x 3 ) . (50) 

The terms in the resulting equation that are purely linear in h a b, contain second derivatives 



of h a b, and do not depend on the background curvature are of the form \d % dih? j — \d % &hi^ 

-6 



i.e. the same combination of derivatives of components^ as appears in ff49l . There are also 
terms which are linear in h a b and contain second derivatives of h a b, which depend on the 
difference between the exact background metric g a °^ and r] ab , and these can be expressed in 
the form U % ^ kl didjhku where [/ y = 0(x 2 ). The terms in the equation which are nonlinear 
in h a b that contain second derivatives of h a b can be expressed in the form diW 1 + Z\ where 
W % converges to zero uniformly on compact sets as A — > and Z\ is uniformly bounded on 
compact sets as A — > 0. The remaining terms in this equation then take the form 87rT 00 + Z 2 , 



6 These terms correspond to the linearization of the scalar curvature of the spatial metric, as would be 
expected from the general form of the "Hamiltonian constraint equation" . 
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where Z2 is uniformly bounded on compact sets as A — > 0. Now multiply this equation by 
fp. Using fH9|) . we see that the resulting equation takes the form 

d%ip{\) = Anf s P T 00 (\) + <W(A) + C(A) + ^U ijkl d i djip k i(\) , (51) 

where u % = \f P W l is of fixed compact support and converges to zero uniformly as A — > 0, 
and ( is of fixed compact support and is uniformly bounded as A — > 0. Thus, (p satisfies a 
Poisson-like equation. Furthermore, the position space version of (|47|) can be written as 

2E V = -— [ dtd^xd^d^ , (52) 

47T J 

which is just the usual formula for (twice) the gravitational potential energy in Newtonian 
gravity! Note that we have not made any Newtonian approximations, nor have we made a 
special choice of "time vector" t a . 

Thus, we see that the resolution of the issue of whether t$ > depends on a competition 
between the positive contribution from the tensor modes and the negative contribution, E v , 
arising from a Newtonian-like gravitational potential energy. We will now show that if 
Too(X) > 0, then, in fact, E v — )■ as A — > and 6 — > 0. Thus, the scalar modes make no 
contribution in this limit, and the tensor modes always "win". 
To prove this, we use ( 15T|) to rewrite E v as 
1 



e <p = it dtd^xipd'diip 

671 



^T 00 (A) + -W(A) + J-C(A) + ^ kl d % d^ kl {\) 

v 47T 47T 671 



(53) 



dt d 3 ^(A) 

The second term can be written as the time integral of 

d^cpdi^ = - J d^xd^J . (54) 
By the Schwartz inequality, we have 



d 3 a;<9i</9u; ! 



< ||<9iy3|| L 2||u/|| L 2 . (55) 



However, ||w*||^2 — > as A — > 0, and we have already noted that ||9j(/9||^2 remains uniformly 
bounded as A — > 0. Therefore, we see that the second term vanishes in the limit as A goes 
to zero. 

To analyze the remaining terms on the right side of (1531) . suppose we could show that 
<^(A) converges uniformly to on compact sets as A — > 0. Then since ((X) is of fixed compact 
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support and is uniformly bounded as A — > 0, it follows immediately that J dt d 3 X(p(\)((\) — > 
as A 0, so the third term in fl53|) vanishes in the limit as A — > 0. On the other hand, 
if (f{\) converges uniformly and T 00 (A) > 0, then the first term is exactly of the form to 
which the above Lemma of this section applies, with A = (p, B = T 00 , and fp being the test 
function with which A(X)B(X) is being smeared. The Lemma then states that the first term 
in ( |53|) vanishes in the limit as A — > 0. Finally, the last term of (|53|) may be re- written as 

dt da; [U m d iV {\)d^ kl {\) + <9^' fc V(A)<9^M(A)] . (56) 
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If converges uniformly to zero, then the second of these two terms vanishes in the limit 
as A — > because U %]kl is independent of A and djipki(^) is uniformly bounded as A — > 
and is of fixed compact support. In contrast to all the others, the first term above does not 
converge to zero as A — > 0. However, it will still vanish once we subsequently take 5 — > 0. 
To see this, use the fact that on the support of f P , U %3kl is bounded by a constant times 5 2 
(see f l50|) ) to write 



dtdxU^d^Wd^kiiX) 



<C5 \\di(p\\ L 2\\djijjki\\L^ < C'5 \\djip k i\\ L 2 ■ (57) 



where the last inequality follows from the fact that ||dj<£>||^ 2 < f H^V^Hx^- However, by f )40|) 
we have 

lim lim ||<9/0fcj|| l 2 = /J? f k i(P) ■ (58) 

Consequently, the right hand side of (157|) vanishes when the limits as A — > and 5—^0 
are taken. Thus, we will have proven that > provided only that we show that <^(A) 
converges uniformly to on compact sets as A — > 0. 

To prove uniform convergence to of y?(A) on compact sets, we note that it follows 
immediately from fH9|) that 

A<p = + X , (59) 



where 



d% X = &&il>ij ■ (6°) 



As already noted above, ipij{\) — and hence ip l i{\) — converges to uniformly as A — > 0. 
Thus, will converge to uniformly on compact sets if and only if xW converges to 
uniformly on compact sets. We will now prove this by "brute force" . 
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The solution to (!60|) is 

Changing integration variables to y = x' — x and integrating by parts, we obtain 

x (t,x) = ~ [ d 3 yd^(t,x + y)^L. (62) 
An J \y\ d 

For any r > 0, we can break up the integral appearing on the right side of this equation 

into an integral over \y\ < r$ and an integral over \y\ > tq. We leave the first integral alone 

but do another integration by parts on the second integral. We thereby obtain 

X (t,x) = ~ [ d'yd^^x + y^ + l- [ dni*^L*(t,x + y)y±) 

+1- f d 3 y^{t,x + y) 6M *- yiyi . (63) 

47r J\y\>r \vr 

Now let -Fi(A) = sup( ta! ) \ipij\ and let F 2 (X) = sup^^ \dkipij\. Then F\ — > and F 2 remains 
bounded as A — > 0. It follows straightforwardly from ( 1631) that 

\ X (t,x)\ <c 1 F 2 (A)r + c 2 F 1 (A) + c 3 F 1 (A)|ln(C/ro)| , (64) 

where ci, c 2 , C3, and C are constants (i.e., independent of A, r , t, and a;). This bound holds 
for all ro and all A. Therefore, as we let A — > 0, we are free to choose ro to vary with A in 
any way that is convenient. Choosing 



r (A) = exp[-l/ v /F 1 (A)], (65) 

we obtain the bound 



\ X {t, x)\<d exp[-l/V^UA)] + C 2 F 1 (X) + C 3 ^/F^(X) , (66) 

from which it follows immediately that x ~ > uniformly as A — > 0, as we desired to show. 
We have thus proven 

Theorem 2. Given a one-parameter family g a b(X) satisfying assumptions (i)-(iv) above, the 
effective stress energy tensor t^} appearing in equation ( fl5l) for the background metric g^ 
satisfies the weak energy condition, i.e., 

t^t a t b > (67) 
for all t a that are timelike with respect to g^J . 
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It should be emphasized that all of the results of this section apply to an arbitrary 
one-parameter family g a &(A) satisfying assumptions (i)-(iv). In particular, no symmetry or 
other assumptions concerning the background metric, g^ b \ were made. However, if FLRW 
symmetry is assumed for gj^ as well as for the (weak limit of) the matter stress-energy 
tensor, then t^j must also have this symmetry. It then follows immediately from 

Theorems 1 and 2 that must have the form of a perfect fluid with P — |p and p > 0. 
In particular, the effective stress-energy tensor arising from nonlinear terms in Einstein's 
equation associated with short-wavelength inhomogeneities cannot produce effects similar 
to those of dark energy. 



III. PERTURBATION THEORY 



In the previous section, we obtained the equation satisfied by the background metric, 
g^ b \ derived key properties of the effective stress-energy tensor t^j , and thereby proved that 
small scale inhomogeneities cannot mimic the effects of dark energy on large scale dynamics. 
However, in cosmology and other contexts, we wish to know not only the dynamical behavior 
of but also the dynamical behavior of the deviations from g^ b \ as this is needed to de- 
scribe the formation and growth of structures in the universe. In particular, we would like to 
obtain the equations satisfied by h ab (X) to sufficient accuracy that h ab (\) can be determined 
to first order in A, i.e., any deviations from an exact solution (over a compact spacetime 
region) go to zero faster than A as A — > 0. As already mentioned in the introduction, if we 
were in the context of ordinary perturbation theory where g a b{X,x) is jointly differentiable 
in A and x, we would define 

A A >° A 



IS 



To derive the equation satisfied by 7^, we differentiate the Einstein equation with respect 
to A, at A = 0. The result is an equation that sets the linearized Einstein operator acting on 
7^ equal to the derivative of the stress energy tensor with respect to A, evaluated at A = 0. 
We would then take h ab {\) = At^. However, in the context of our framework, g ab (X,x) is 
not differentiable in A at A = 0, so we cannot even define a notion of a "metric perturbation" 
by differentiating g ab (X,x). 

Of course, the (exact) equation satisfied by h ab (X) is simply the equation obtained by 
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substituting g a b(A) = g^b + h ab {\) into Einstein's equation (fl2~|) . However, this is not any 

more useful in practice than simply asserting that g ab {X) must be a solution of Einstein's 

equation for all A; if we could solve Einstein's equation exactly, there would be no need to 

develop a perturbative formalism. The key idea needed to obtain a more useful version of 

Einstein's equation is that, although nonlinearities may be important on small scales, there 

should be a simpler, linear description on large scales. The key idea needed to implement 

this description is the observation that although the ordinary (pointwise or uniform) limit of 

[gabW ~9 < fb]l^ does n °t exist in the context of our framework — g a bW is not different iable— 

there is no reason why the weak limit of this quantity cannot exist. 

Thus, a natural generalization of the conventional linearized metric perturbation is 

(L) ,. flw>(A) - g^b / Rf y. 

Tab = w-lim . (69) 

A^O A 

Here we have replaced the ordinary limit of ( 168]) with a weak limit, which "averages away" 
the small scale inhomogeneities. This quantity thereby corresponds closely to the notion of 
the "long wavelength part" of the metric perturbation that appears in other analyses (see, 



e.g. 



28]). We will discuss this further in section V below. The remainder of the perturbation 
will be denoted 

fcS?(A)sMA)-A7£\ ( 7 °) 
and will be referred to as the "short wavelength part" of the deviation!!! of the metric from 

gf b \ In subsection IVB, we will argue that, under suitable Newtonian assumptions in cos- 

( S) 

mology, to leading order in A, h ah depends only locally on the matter distribution and is 
well approximated by a Newtonian gravity solution. We emphasize that, within our frame- 
work, short and long wavelength perturbations have a very different character. The long 
wavelength part of a perturbation has a well-defined description in the A — > limit, namely 
7oi . On the other hand, the short wavelength part, h^(X), is defined only for A > and 
has no description in terms of a limit as A — > 0. 

We can obtain an equation for by taking the difference of the exact Einstein equation 
( IT2|) for (? a fe(A) and the background Einstein equation ( FT5|) for g^ b ' , dividing by A, and taking 
the weak limit as A — > 0. However, in order to ensure that 7^ is well defined and satisfies 
7 If we were to consider higher order perturbation theory, then we would also subtract from h a t, higher 

order in A "long wavelength" contributions to define h^(X). However, we shall only be concerned with 

first order perturbation theory in this paper. 
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a well-defined equation, we must append to assumptions (i)-(iv) of section II the following 
additional assumptions on our one-parameter family 5W>(A): 

(v) There exist smooth tensor fields 7$, ^ cdef , v^Le and ^abldefgh 011 M such that 



(a) 



(b) 



(c) 



(d) 



W A -hmi/aA)=7^, (71) 



w-lim — 

A^O A 



V (a ^I3(A)V 6) ^ ; (A) - to/ = ^Lf , (72) 



w-lim — 

A^O A 



h ( b S c \\)V a h[ s J(\)]=uZ e , (73) 



w-lim — 

A^O A 



h%\\)V a hlf(\)V b h%\\) = ^> cdefgh . (74) 



(i) 



In the following, we shall assume that our one-parameter family g a b{^) satisfies assumptions 
(i)-(iv) of section II together with assumption (v) above. In this section we will make no 
additional assumptions about g ab (X), so all of the results obtained in this section should hold, 
e.g., for self-gravitating gravitational radiation in a background without any symmetries. In 
section IV, we shall specialize to the case of Newtonian-like cosmological perturbations off of 
a background metric with FLRW symmetry, and will make numerous additional assumptions 
and simplifications. 

The newly defined first order backreaction tensors, H a l bcde f, v< ab\de an d ^abidefgh^ possess 
certain tensor symmetries as a direct consequence of their definitions. Clearly, as with the 
zeroth order quantity fi abcdef , we have $ b){cd){ef} = ^ bcdef and ^ cdef = ^ ba ] efcd . However, 
in contrast with fi a bcdef, the symmetry under interchange of the first two indices had to be 
built directly into the definition of A^de/' ra ther than derived. Indeed, 

w-lim i [V a h cd (\)V b h ef (\) - fi a bcdef] 

A— >0 A 

= Vabcdef + w-lim ^V [o /l H (A)V 6 ]/le/(A) 

= Vabcdef + \ [ V I« (Vl(A)V b ]MA)) - MA)V [a V 6 ]/l e /(A)] 

= Vabcdef + ^[aVqldef ~ 2 \ k cdW R ab( J 'h f )g (X) 

= Vabcdef + ^[aV { b]ldef ■ ( 75 ) 
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What this calculation also illustrates is that, because of the factor of 1/A, we may no longer 
freely drop total derivative terms when taking weak limits, and in general we will pick up 
terms of the form V a ^ e/ . It follows also that ^)(«te) = "abide and "abode = ~ v SLc- Finally, 

U ab\cd)(ef)(gh) = U abode fgh^ W abode fgh = ^bacdghef alld 

Lo {1) -Lo {1) =ui (1) (76) 

abode f gh bacdefgh baefcdgh abefcdgh ' V / 

We also note that definitions (b) and (c) would be unchanged if h ab were replaced by h a b, 
and (d) would be unchanged if any h ab which is being differentiated were replaced by h a b- 
We now subtract the background Einstein equation (fT5|) from the exact Einstein equation 
(TT2]) . divide by A, and then take the weak limit as A — > 0. A very lengthy calculation, 



performed with the help of the xAct tensor manipulation package 30] for Mathematica, 
yields 

v c v (a <) - \ v<v c7 ^ - \ v a v 67 ^ c c - l g $ (W 7 £> - v c v c7 ^ d d ) 



1.1 1 



6 6 g | 6 6 616 

2^abc de ~ 2^c(ab)d e ~ 2^ ( ab )ecd ' ~^ i cdab e — Tj^^cda be ~ A*c abde + A*c e(a&)d 



+ eabcd — 7^ eacM g^a^ { _ ^cd If ~ ^cd ^ ef + d^e/ — ^/i eC( /j + 2/i ^ ce( ^|^ 



-T (1) 4- -n (0) i-i/O de _ e , 9 (l)cd e l,I (l)cd 

J a6 + gi/ab IP c de A* c d e + Z P c de] ^ acbd 

1 ..(l)c d , 1(1) cd _I,.(1) C d +-,/( 1 ) c d - l/yW^ 

6d ^P ab cd (ab)c d ' ^ cab d abcd 



+ {^ 1)C c de dJ } + ^cVef + ^jjf + ^^ef ~ ^cYe, ~ deof) 



1(1) cd e , 1(1) c d e 1(1) ede _ 1 (l) c d e _ 1 (l) c de _ 3 (1)c de 

2^ ab c de ^ 2 (a |c| 6)d e 2 ^ a ^ cde 8 ^ d 6 8 ^ 4 C 



, I (l)c de , 1 (l)cd e _ ^ (l)cd e , 1 (l)cd e , 1 (l)cd e _ 1 (l)cd e 

2 c adbe 4 f^^d ce 4 cdab e O cc ' a ^ e 2 C a ^ e 2 C ac ^ e 

e _ e , 1(0) /oy .,,(l)c de 4. V /V^ d e l - Iv, 7/ (1)C d 

+ 2 rf abce 2 ^ aebC 4 t ZVdZ ^ c e+ V e V cdj 4 V ( a ^ 6)c d 

4-iy i/ 1 ) cd _ Iv 7 y( 1 ) c d - X/jiA^ cd + V^W c d -u Iv^yW d (77^1 
+ 4 c ("b) 2 ab d (afe) c d afcc 2 c(af>) ' 

Here, we have written 

T a « ^ w-lim Tab{X) " ^ . (78) 

This weak limit exists by virtue of assumption (v) and the fact that g a b{X) satisfies Einstein's 
equation. 

27 



Our equation ( 177|) for 7^ takes the form of a modified linearized Einstein equation. The 
terms on the left side are linear in 7^ and, in addition to the usual terms appearing in the 
linearized Einstein tensor, contain terms proportional to Habcdef- The right side contains, in 
addition to the matter source , numerous "effective source terms" arising from y^bcdefi 

1 'abode anC ^ U abode fgh' 

Some relations between jtt£2xfe/' v abode anc ^ u abode fgh can ^ e derived from Einstein's equa- 
tion. We previously derived the relation a aeb e cd = (see (I2B1 ) by starting with Einstein's 
equation in "Ricci form" ffl9|) . multiplying it by h e f(X) and taking the weak limit as A — > 0. 
In a similar manner, if we multiply ( fT9l) by /i^p (X)hfj (A) / A, we obtain the following equation 
satisfied by d£L*fgh- 

mWs 9 

acdefbg 1 bcdefag gcdefab abcdef g 

^ W aefcdbg ' ^ befcdag M gefcdab w abe/cd g ~~ u ' I ' ^ ) 

In deriving this, we used the fact that 

w-limi/»2 ) (A)/»S ) (A)T e/ (A) = ) (80) 

A— >U A 

which follows from our lemma from section 2 and the assumption that T a b(A) satisfies the 
weak energy condition, in the same way that we showed n a bcd = (see (1271)). 

Similarly, if we subtract the background Einstein equation in "Ricci form" from (119p . 
multiply by h ( cd \X)/X, and take the A — > weak limit, we can derive an analog of (128"]) 
satisfied by ^ cdef , 



I, /(l) e -L _ „(l)e , , 

2^ abed e T 2" caiab r 1 (a\cd\b)e 



= 87r/t (1) cda6 - 47T^° ) K (1) c / e + ^7 (L)e/ (/^abe/cd + Vefabed ~ ^e{ab)fcd) + ^ W abed^ ef 

^2 afe cde f 2 ^ a ' Cd '^ e / (° l e l b )f cd 4 ecdafe / 2 6Cda 7 2 cda/be 

+ 2^ efabed+ 2 i a 6)aie+2 (" b)ecd + g e ( a fe)cd 2 Ve abcd ' ^ ^ 

In this equation we have defined the quantity 

4L = w-lim^i?(A)T cd (A). (82) 

The existence of this limit is guaranteed by Einstein's equation together with our other 
assumptions. Note that since h ( f b ) (X)/X does not converge uniformly to zero as A — > 
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(although it is uniformly bounded), the weak energy condition does not imply the vanishing 

° f K abcd aS Jt did for K ^cd- 

We can simplify the above equations as follows. As with the background case, define 
a abcde f = $\[ab]\d]ef and PjLef = ^abcd)efi and express the equations in terms of these 
quantities. A similar breakup is also possible for ^bcdefgh- First, split it into the parts 

(i) 

(ab)cdefgh 

l-<-tr II V Tim hi mri-i 

abcdefgh [ab]cdefgh' 



with are symmetric and antisymmetric in the first two indices, 0JabfJ e f gh = ^fa\)cdefgh and 

UabSefah = "UMcdefah- B ^ MM , We haVe 



U 'abcdefgh ~ U abefcdgh ' \° 6 ) 



The symmetric part is then decomposed into u) < ab ' c J e j gh — ^\ e \[ a \ C( i\b}\f]gh aIld ^abflfgh 
u (ab\cd\ef)gh' w ^ n inverse transformation 



(1,S) _ _4 / (l, a ) _ 

^abcdefgh ~ 3 ^eHcdlfO/g/i a '(/(a|c<i|fe)he/ ^(elcdl^/iab 

+ a; (1>/3) -£j (1,|S) 

'abcdefgh ' abcdghef efcdghab ' 



(84) 



Substituting for our new quantities, (J79|) can be rewritten as 



^'^aW/ + ^Jefbgcd = , (85) 



whereas (IHTl) becomes 



„(1) e _ 4 (1) _ 27rf7 (°) K (l) e 4. „ . , 4. 1 W (M) e / _ jL(M) e / 

« a becri — 47rK cda6 ^9ab K cde + a aebfcdl + . W (a b)ecd f n U tab). 



L,(M) e / L,(l,«) e f (i, a ) e / (l, a ) e / _ 1 (l,«)e/ 

^ aeb/cd ^2 a cdb ef ' u (a |cde| 6)/ ' W a e b/dc ^ feae/ric 



(a 6)ec<2 / 2 ( a b ) cde f 

f 4 w (l.«) e / _ jL(L«) e / 
cde| b)f^ U a e bfdc ^ I 

+ -^^ 1 \ ) cd e e + \^^ 1)e b)ecd + \Ve^\ ab) \ d - \v^\ bcd . (86) 

Finally we can use ( 185]) and (IHBI) to simplify our version of the linearized Einstein equation, 
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([77]) . resulting in 

v c v (a7 £ - iv<v c7 i? - ±v.v t7 (L)B fl - (V v^ d L) - v c v c7 ^ d 

= 8vtT q « + a^ a \ d cd - 2™« ab c c - 8™« ( / 6)c + 2^<® {*« c /d - /c^} 

_1 (1,A) c d e 1 (M) c d e 1 (l,A)cd e , 1 .(l,A)cd e 

^ W (a b) cd e + 2 ( a b ) c de 8 ^ de 4 6 

+ — n {0) lu,^ cd e f + 2tJ l ' A ^ cd e f \ - ced - uM c de - (J 1 '") 

16 I C de f cdefj ZW (a b)dce U (a b)d ce U (a b)cde 

, 1 (l,a)cd e , 1 (l, a )od e , 1(0) f (l,a)c edf _ o, j; (l,«) c e/d 1 
oid e + ^ W bacde^ g3ab \ U dec f ZU} dec fj 

_Iy, t/( 1 ) c d + iy. z/ (1)c d +ivi/ (1)c d - -V^ (1) Cd + -V,v {1)c d 

2 ( a l)H + 4 V K fe)cd+ 4 V c I/ ab d 4 Vd ^ (06) c + 4 Vc ^ a6c 

^ ,/( 1 ) c d -L 1X7 ,/W c de -L Y7 ,/( 1 ) c d e l (R7\ 



ce 



Equations (J8"5j) . (|H5|) and f l8"7j) describe the long wavelength perturbations. It should be 
noted that, just as (3 a bcdef was absent from our background equations, f3 a bcdef, Pabcdefi anc ^ 
u ab'cdef g h are an absent from our perturbation equations. 

Equations (1851) . (I8"6"j) and (18 7p have been written down in an arbitrary gauge. We will make 
a specific choice of gauge in subsection IVA below, but for now we note that we can apply 
any one-parameter family of diffeomorphisms, <ft\, to 5w,(A) that preserves conditions (i)- 
(v). Burnett 26J has analyzed the properties of gauge transformations associated with one- 



parameter families of diffeomorphisms that are not smooth in A. Here, we simply note that 
any smooth, one-parameter group of diffeomorphisms cf)\ generates gauge transformations 
that are easily seen to preserve conditions (i)-(v). Under such gauge transformations, it is 
not difficult to see that 7 ^ — > 7^+^^ , where £ a is the vector field that generates <p\ and 
C denotes the Lie derivative. Thus, 7( ^ h as the same gauge freedom arising from smooth 
0a as in ordinary linearized perturbation theory. This freedom can be used to impose the 
same types of gauge conditions on j^' as in ordinary linearized perturbation theory. It 
is also not difficult to see that h a S h (A) — > <^/t^/ (A) + j a b(A), where j a fe(A) = 0(X 2 ) and is 
jointly smooth in A and the spacetime point. By using this gauge transformation property 
of Kb W> it; is Possible to show that fi a bcdef, v^,, ^abldefgh^ and K abcd are g au S e invariant 
under gauge transformations arising from smooth (J)\, whereas A*SLfe/ — rabcdef + £>zHabcdef 

We turn our attention now to the short wavelength perturbations. Without making any 
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approximations it is straightforward to write down an equation satisfied by h^\X): Simply 
substitute g a b(X) = + Xry^' + nf h \X) into the exact Einstein equation. We may write 
this equation in the form 

G^V 0) , ^ s \X)) + AA<?(A) = 8nT ab (X) ~ GM 0) ) ~ A ff £ } - A^V 0) ,7 (i) ) 

oo 

-AAT^-E^^^T^ + ^CA)), (88) 

n=2 

where we have grouped linear terms in h^(X) on the left hand side. Here, G^(g(°\ X'y^ + 
h^(X)) denotes the nth order Einstein tensor expanded about g^J of the perturbation 

Unfortunately, it does not appear possible to simplify (IBB]) to obtain suitable approximate 
solutions without introducing additional assumptions. Of course, if we do not simplify 
(|88|) . then we have not made any progress beyond asserting that we must solve Einstein's 
equation. In the next section, we will introduce additional assumptions relevant to the case 
of cosmological perturbations and argue that to obtain an accurate description of the metric 
to order A, we may replace (188]) by the equations of Newtonian gravity with local matter 
sources. 

For the remainder of this section, we shall compare our general analysis to that given 
by Isaacson 23|, |24J and others (see, e.g., [25]), who were interested in describing the self- 
gravitating effects of gravitational radiation. We therefore restrict attention to the vacuum 
case (T ab (X) = 0). These authors work with the quantity h ab (X) rather than introducing 
Habcdef- Suppose one is merely interested in determining the background metric g^ b \ i.e., 
one is not interested in obtaining an accurate (to order A) description of the deviation of 
the metric from g^ b \ Then one would need only to calculate h ab (X) to sufficient accuracy 
that one could determine fi a bcdef and, thereby, t® b (see ( IT6]) ). In particular, one would not 
be interested in computing "f ab \ so one could ignore the equations we have derived above 
for 7^. Furthermore, in order to obtain h^f^X) to sufficient accuracy, it appears plausible 
that one could make 0(1) modifications to flHHl) as A — > and still determine h^\X) to 
sufficient accuracy, provided that these 0(1) modifications have vanishing weak limit. Here, 
by the phrase "it appears plausible" we mean that we believe it is likely that one could 
introduce additional reasonable assumptions on the one-parameter family g a b{X) so that 
these modifications to (188]) could be made without affecting g^ b \ The reason for this belief 
is that 0(1) error terms in (IBB]) should — under suitable further assumptions similar to ones 
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indicated at the end of subsection IVB below — give rise to 0(X 2 ) errors in h^\\), which 
should not affect Habcdef- 

A candidate modification of flHHj) in the vacuum case would be to drop the entire right 
side of this equation (together with the term Ah ab (A) on the left side), to obtain simply the 
linearized Einstein equation for off of g^ h , 

G , il ) (^ (0 U (S) (A)) = 0. (89) 

However, if G a b(g^) 7^ 0, the linearized Einstein equation off of g^ b does not appear to have 
an initial value formulation, so this modification of ( 188]) is probably not suitable. (Note that 
the linearized Einstein equation off of a non-solution is not gauge invariant, so one cannot 
employ a choice of gauge to simplify the equation and/or put it in hyperbolic form.) A 
better candidate modification would be an equation of the same form that the linearized 
Einstein equation would take in the Lorenz gauge if perturbed off of a vacuum spacetime. 
Since the linearized Einstein equation off of a non- vacuum spacetime is inconsistent with the 
Lorenz gauge condition, there is no unique choice of such an equation — in particular, one 



can add new terms involving the background Ricci tensor — and, indeed, Isaacson 23j, [24 ] 



and Misner, Thorne, and Wheeler 



25] give slightly different forms of the proposed equation 
(compare Eq. (5.12) of |23j with Eq. (35.68) of j^). In fact, since terms involving the 
product of the background curvature with h^ b are 0(A), it would appear simplest to drop 
all of these terms and work with the wave equation 

V C V C /^ = , (90) 



where h!f b = fif b — \g^hS s ^ c c . As with the equations used in 23] and 25], this equation 
is inconsistent with the Lorenz gauge condition V a h ab = 0. However, if one constrains the 



initial data for h ab for solutions to fl90|) so that V a h ab and its first time derivative vanish 
initially, then the Lorenz gauge condition should hold to 0(A) at later times (in compact 
regions of spacetime). Thus, the solutions to fl90|) with these initial data restrictions should 
satisfy Einstein's equation ( I8"8"j) to the desired 0(1) accuracy. It should also be possible to 
impose the gauge condition h^ a a = to 0(A) accuracy. 

If the above arguments are correct, then in order to obtain the possible background 
metrics g^ b \ it should suffice to simultaneously solve ( !T5|) and ( 190]) (or equivalently, Eq. (5.12) 
of 23] or Eq. (35.68) of |25[) with appropriate restrictions on initial data, to obtain g^J and a 
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one-parameter family h ab (\) satisfying our conditions (ii)-(iv), where, in (j!5p . t ab is given by 
( IT6|) and Habcdef is given by (JHJ) with h ab (\) replaced by our one-parameter family of solutions 
to f )90|) . One could then attempt to establish properties of t^ b — in particular, the vanishing 
of its trace and the positivity of its energy density — by working with the expressions for 
it in terms of h ab (A) in a particular gauge. Needless to say, numerous extremely murky 
mathematical issues arise if one proceeds in this manner. Our analysis of section II, as well 
as the work of Burnett |26| in the vacuum case, proved rigorous results about tf b without 
introducing any approximate equations satisfied by (A), and thus completely bypassed 
these murky issues. 

Finally, we note that if one wished to know the deviation of the metric from g^ b to 
0(A) — as would not necessarily be of interest in studying gravitational radiation but is of 
considerable interest in cosmology — then, of course, it would be necessary to know 7^. 
Although we argued above that one has considerable freedom in modifying the equation 
satisfied by h^ b without affecting g ® b , it is clear that one cannot make any modification to 
the equation ( 1BT1) satisfied by 7^ without introducing 0(1) errors in 7^ and thus 0(A) 
errors in g a b(ty- Furthermore, in order to calculate the effective source term (y^ a c b d cd in 
fl87|) . it appears that one would need to know hf b to 0(A 2 ), in which case one could not 
drop the quadratic terms in h ab in (I88J) . Thus, in the case of self-gravitating gravitational 
radiation, if one wished to know the deviation of the metric from g^ b to 0(A), one would 
have to solve ( IT5l) . ( 1H71) . and some suitable simplification of ( )88l) . This would comprise an 
extremely complicated system. Fortunately, in the case of cosmology, we will now argue 
that, under additional assumptions, significant simplifications occur. 



IV. COSMOLOGICAL PERTURBATION THEORY 



Up to this point we have not assumed any symmetries or other special properties of 
the background metric g^J . We also have not made any restrictions on the matter con- 
tent, T a h(A), other than that it satisfy the weak energy condition, nor have we imposed 
any restrictions on the perturbations. In this section, we will be concerned with the case 
of main interest in cosmology, where g^ b has FLRW symmetry, there is negligible gravita- 
tional radiation content (in particular, t^ b = 0), and the matter content satisfies suitable 
Newtonian assumptions. We will argue that, under these assumptions, to leading order in 
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A, hJl (A) is given by local Newtonian gravity, i.e., in a neighborhood of any point x, hj^ 
can be calculated to sufficient accuracy using Newtonian gravity, taking into account only 
the matter distribution within a suitable neighborhood of x. (The effects of more distant 
matter are taken into account by 7^.) With our "no gravitational radiation" assumption 
for the background and our Newtonian approximation for , our equation ( 1H71) for 7^ 
simplifies considerably, yielding the linearized Einstein equation with an additional effective 
source that agrees with recent results of [28( (see also [2l|). 

In our analysis, it will be important to make a convenient choice of gauge. Since we are 
taking nonlinear effects at small scales into account, we cannot simply impose the usual 
cosmological gauge choices for perturbations, i.e., we must make our gauge choice at the 
nonlinear level. When studying non-linear perturbations off of a flat background it is often 
convenient to work with "wave map coordinates" (usually called "harmonic coordinates" 
in the literature), particularly in the context of the post-Newtonian expansion [jy]]. Since 
our background metric is not flat, the usual definition of wave-map coordinates is not 

impose a generalized wave-map gauge condition with respect 



convenient, but we may insteac 
to the background metric g^} 



32]. 



Our gauge choice will be introduced in subsection A in the context of a general background 

metric g^J (i.e., without assuming FLRW symmetry). In subsection B, we restrict to a 

( S) 

FLRW background, we make our Newtonian assumptions, and argue that h ab is given by 
local Newtonian gravity. It should be emphasized that the arguments of section B have the 
character of plausibility arguments rather than proofs. Finally, the simplifications to the 
equation for 7^ will be obtained in subsection C. 



A. Generalized wave map (harmonic) gauge 

Given a one-parameter family of metrics g a b(A) on our spacetime manifold M that satisfies 
our assumptions (i)-(v), we may apply any one-parameter family of diffeomorphisms, 0(A) : 
M — y M, which preserve these conditions, where without loss of generality, we may assume 
that 0(0) is the identity map. As already noted in the paragraph below ( 1H71) . it is clear 
that any 0(A) that is jointly smooth in A and the spacetime point x will preserve conditions 
(i)-(v), but there also should exist a wide class of 0(A) that are not smooth in A that 
preserve these conditions. The properties of such 0(A) and the transformations that they 
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induce on Habcdef were analyzed by [26( (under his assumptions, which differ slightly from 
our assumptions (i)-(iv)). Unfortunately, although Burnett's analysis can be used to prove 
important properties of gauge transformations, such as the invariance of t^j under all allowed 
gauge transformations, it is very difficult to prove any existence results that establish that 
specific gauge conditions can be imposed on an arbitrary one parameter family of metrics 
^ a b(A) satisfying our conditions. 

In this section, we will assume that we can impose the "generalized wave-map gauge 
condition" on the metric <7 a &(A), namely 

g ab (\)C c ab (\) = 0, (91) 
where C c ab {\) is given by (flQj) . Note that this condition depends upon the background 
metric, g^ b \ since the derivative operator, V a , of gf b appears in the definition of C c ab (X). We 
can impose the gauge condition ( 19~TT) on g a b(\) by applying the diffeomorphism x M — > </> M (A, x) 
to g a &(A), where x M are arbitrarily chosen coordinates on M and M satisfied 
1 d 

This equation is a (nonlinear) wave equation, so we can always find (local) solutions. To 
see that solutions to this equation give rise to the condition ( )9T]) . we note that if we use 
M (A,x) as coordinates for the Ath spacetime, then we may replace d(j) a /dx^ by 5°^, and 
( |92l) reduces to 

=9* (V^iWg^w) + r<°> v*r(A) 



= r(A)rv(A)-^(A)r(°)v 

= -g^{\)C\ v . (93) 

We will refer to the coordinates M (A,x) as generalized wave map coordinates for g a bW 
relative to the coordinates x M for gfj. Note that in the case where gfj = rj a b and x^ are 
Minkowski coordinates, the generalized wave map gauge condition reduces to the condi- 
tion g^ u (X)T a ^ u (\) = 0, and the coordinates M (A,x) satisfy a linear wave equation. Such 
coordinates are usually referred to as "harmonic coordinates" in the literature. 



The diffeomorphisms defined by (|92|) do not depend on the choice of coordinates £ M , since this equation 
can be derived from the coordinate invariant action 
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Although we can always (locally) solve ( 192|) and thus (locally) put each g a b(X) in our 
one-parameter family in wave map gauge, we have no guarantee that the resulting new one- 
parameter family of metrics will satisfy our conditions (i)-(v). In the following, we shall 
simply assume that this is the case, i.e., that we have a one-parameter family of metrics 
g a bW that satisfies conditions (i)-(v) as well as our gauge condition ( 19~T|) . This corresponds 
to a strengthening of our assumptions. 

When the generalized wave map gauge condition is satisfied, it is very convenient to work 
with the variable 

r (A) = g^ ab - y^^(A) • (94) 

instead of h a b(X) = g a b{\) — g% ■ (Note that \J g(X) / g^ is the proportionality factor between 
the volume elements of g a ft(A) and , so this quantity does not depend on any choice of 
coordinates x M on M.) We have 



W 6 (A) = - J^v^(A) - J^(x)cr A 



g(0) 



where on the first line we used the fact that 



(v b g ab (\) + \g a \\)g cd {\)V b g cd {\) 
V c (A)C" 6c 



3(A) 



(95) 



**=5V.bgg$). (96) 

Note that in linearized gravity, h afe reduces to the trace- reversed metric perturbation, i.e., 
h a6 (A) -> h ab (\) = h ab (\) - |g (0)ab /i c c (A), and (jHSD reduces to the Lorenz gauge condition 
V ' a h ab = 0. One should keep in mind, though, that beyond lowest order in A, h a6 (A) is not 
the trace-reversed metric perturbation. 

We now express the exact Einstein equation in wave-map gauge in terms of f) ab and the 
background derivative operator. Starting with (j!2p . we use the background equation ffT5l) . 
the gauge condition Vbf) a6 = 0, as well as the fact that 



V a g bc (\) = (v a ^ c - l -g bc (\)g de (\)VX e ) , (97) 
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to show (after a long computation) that 

V c V c f) afc - 2R a cd \g (0) )t) cd + 2R} a {g {Q) )\) b)c - R(g (0) )t) ab - g [0)ab R c d{g {()) W d 
-2G^(°))^V + 2A (V - ^ (0) ^2V 

= -16vr^ (T ab (X) - T (0)afe + i ab (X) - t {0)ab ) . (98) 
Here, the terms that are non-linear in \) ab have been absorbed into the quantity, 
167r^t a6 (A) = -2G ab (l - 9 j§- - g<$ffA + iU<? (0) )f) c V - 2Rj a (g^)^ d 
-i) cd V c V d i) ab + V d t) ca V c t) db + g^{\)g cd {\)V e ^ a V f \) db 
+ l -g ab {X)g cd {X)VJ) fc Vft ed - 2g cd {X)9 f[a {^e^ b)c V^ ed 
+i {2g a \X)g b \\) - g ah {\)g 9h {\)) {2g cd {\)g ef {\)-g ed {\)g cf {X))V h ^ d V g ^ 

-2A fi^- [g ah {X) - g^ ab ) + \) ab - ^ (Q) ^£V) • (99) 
Equation (195)) is of the form 

j^ ab (f)) = -lQnS ab . (100) 

where J2f a6 takes the form of a linear wave operator acting on \) ab . Consequently, re- 
introducing our (arbitrarily chosen) coordinates x^ of the background spacetime, we may 
rewrite (I98p in the following equivalent integral form: 

\) a \Xx) = 4 f G vct a \, u , (x, W (A, x') J-gW (x') d V + f,S£a(A,x) , (101) 

JM 

where G Iet „' u > (x, x') is the retarded Green's function for Jzf ab and f)^ m is a solution to 
~^ ab (l)hom) = 0. We emphasize that (11011) is not a solution to ( 1981) since the source S 1 "^ 
depends on f)°^. Rather, it is simply a re-writing of ( 1981) in an integral form. 



B. Local Newtonian gravity 

For the remainder of this section, we restrict attention to the case where the background 
spacetime g^ has FLRW symmetry. It will be useful to work in coordinates where the 
metric components are nonsingular. Thus, instead of the more common choice of polar-type 
coordinates, we will write the background metric in the form 

ds^ 2 = -dr 2 + a 2 (r) (l + k(x 2 + y 2 + z 2 )/A) ~ 2 [dx 2 + dy 2 + dz 2 } , (102) 
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where k = 0, ±1, depending on the spatial curvature. 

The "long wavelength part" of the leading order in A part of f) a6 (A) is given by 

7^° 6 = w-lim ^ = 7 < i > a6 - V )o6 7 (L)c c . (103) 
It follows directly from fl95|) that ^ L "> ab satisfies the Lorenz gauge condition 

V a 7 (L)afe = . (104) 
The "short wavelength part" of f) a6 (A) is given by 

^ ab (\) = ff 6 (A) - A 7 (L)ab . (105) 
By O, it satisfies (in the notation of (jlDDjl ) 



S ab (X) - A w-lim ^ ' 



(106) 



where the weak limit appearing in this equation exists by virtue of assumption (v) of section 
III. Note that S ab (X) still depends on the full perturbation, i.e., f) ab cannot be replaced by 
expression for S ab . Equation f TT06l) can be rewritten in integral form as 

tf 3 )»f>(\, x ) =4^ G vct a ^, u ,{x,x') (V v (A,x') -Awjm ^^ j \]-9<®{x') dV 

+C!n Q/3 (A,^) ■ (107) 

Our aim for the remainder of this subsection is to argue that — in the absence of grav- 
itational radiation and under suitable assumptions concerning the behavior of the matter 
distribution T ab (\) — to leading order in A, i)^ ab (X) near point x is described by Newtonian 
gravity, taking into account only the matter distribution in a suitable local neighborhood of 
x. However, in order to derive this conclusion, we must make significant additional assump- 
tions about our one-parameter family g a b(X), and severe difficulties arise if one attempts 
to formulate these assumptions in a mathematically precise manner. The reason is that, 
although there are simple, precise limits that one can take of general relativistic spacetimes 



to obtain Newtonian gravity |33M35| . these limits would not be compatible with our assump- 
tions (unless h a b = 0(A 2 )), and thus would not be suitable for our use. We believe that it 
should be possible to concoct a mathematically consistent set of assumptions that would 
enable us to rigorously justify our conclusions below, but we do not see a simple and/or 
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elegant way of formulating such assumptions, and we do not feel that it would be illuminat- 
ing to attempt to derive our results from a complicated list of technical assumptions whose 
intrinsic plausibility is not much greater than that of the conclusions we wish to draw. Thus, 
in the discussion below in this subsection, although we will clearly indicate the nature of 
the assumptions that are needed, we will not attempt to formulate all of our assumptions 
in a mathematically precise manner, and we will thereby resort to plausibility arguments to 
obtain our conclusions. 

We are interested in obtaining t)( s ' ab (X) near a point x at a time roughly correspond- 
ing to the present time in the actual universe. We first argue that although the retarded 
Green's function integral extends all the way back to the "big bang", it should suffice to 
integrate only over the "recent universe" (corresponding, say, to z < 1000 in the present, 
actual universe). There are two reasons why this should be so: (1) The universe is ex- 
pected to be very nearly homogeneous and isotropic in the distant past, so the source term 
(S ab — A w-liniA'-j-o [S ab / A'] ) should be negligibly small. (2) The nature of the retarded 
Green's function in an expanding universe is such as to make the influence of distant sources 
small (on account of redshift and intensity diminution). Similarly, we assume that the grav- 
itational radiation content of the present universe arising from the "big bang" is negligible. 
Consequently, we discard the last term, i)^J m a ^, in (j!07p . Thus, our integral relation for 
§(s) ab becomes 

f)^(A,x) =4 J G^^x,**) (s^\\,x')-\w-\m, S ^ { ^ X '^ f 3 (»)(x')dV, 

(108) 

where W is the (compact) region corresponding to the "recent universe" (i.e., z < 1000 in 
the present, actual universe). 

Next, we argue that, in order to calculate f^ 5 )"^ 3 to O(A) at x, it suffices to perform the 
integral in ( I108P only over a small neighborhood V of x. In other words, we argue that — for 
any fixed neighborhood, V, of x — as A — > 0, the contribution to the integral in d 1 8 [) from 
the region x' G W \ V should vanish faster than A as A — > 0, i.e., 

J wv ( ^f'^ - w 7 lim gMV W g ') ) yC^^dV . (109) 

To see this, we note that if G T J*^ , v ,{x, x 1 ) were smooth (and if the sharp boundaries of 
the region of integration were replaced by smooth cut-off functions), then we would be 
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integrating the quantity (S^/X — w-limy^olS^ / \']) with a test function. Since, clearly, the 
weak limit as A — > of this quantity is 0, it follows immediately that (I109p would hold. Of 
course, G Tet , v ,(x,x') is not smooth, so (I109p cannot be expected to hold without further 
restrictions on S ab . However, in a normal neighborhood of x, the retarded Green's function 
G ret u'v'( x ) x> ) f° r the linear operator Jzf ab takes the form 

GJ^^M = U a ^Ax,x')5 + (a) + V a %,(x,x>)6 + (-*) , (HO) 

where U al3 ^ u ,(x,x') and V a/3 ^ u ,(x, x') are smooth bitensors and a is the squared geodesic 
distance between x and x' (see 36( for details). Thus, apart from the singularity at x' = x 
(which is excluded from the region of integration in (j!09p ). the singularities of G ret a ^ , v , (x, x') 
are of the form of a restriction to the past lightcone (i.e., S + (a)) and a cutoff at the past 
lightcone (i.e., 6+ (—a)). If S ab is not rapidly varying with time (as should be the case under 
our Newtonian assumptions below), these singularities should be quite benign, so it seems 
not unreasonable that (I109P will hold under suitable assumptions. 

We have argued that (I109p should hold for an arbitrary neighborhood V of x. However, 
the vanishing of the contribution to f)^ a/3 (A, x) from outside of V to 0(A) holds only in the 
limit as A — > 0, and the smaller we take V, the smaller we must take A in order to get a 
good approximation to i)^ a/3 (\, x) by integrating only over V. At any finite A, we cannot 
take V to be arbitrarily small and still get a good approximation to f)^ a ^(A, x). How large 
must we take V at finite A? 

To propose an answer to this question, we restrict consideration to the case where the 
matter content satisfies suitable Newtonian behavior (at least in the "recent universe"). We 
shall assume that as A -> we have T 00 (A) = 0(1/A), T 0i (A) = 0(1/A 1/2 ) and T^(A) = 0(1), 
so that, for small A, the energy density is much greater than the momentum density, and 
the momentum density is much greater than the stress. In addition we shall assume that 
spatial differentiation of components of the stress-tensor results in blow-up as A — > that is 
a factor of A" 1 faster than the undifferentiated components (so, e.g., diT 00 (X) = 0(1/A 2 )), 
but that time differentiation results in a blow-up of only a factor of A" 1 / 2 faster (so, e.g., 
doT o(A) = 0(1/A 3 / 2 )). 

We now introduce the notion of the scale of homogeneity, Z(A,r ), at cosmic time r as 
follows. First, we define a fiducial window function which is to be used for averaging. Fix 
a time interval At t and let Vn tXQ to be the "cylinder", centered at Xo, of "height" At 
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and proper spatial radius R. Let Xr,x (x) be a smooth non- negative function which is equal 
to one on Vr )Xo , and falls rapidly to zero outside of this region. Let H denote the "Hubble 
volume" relative to x at time r , i.e., the ball of proper radius equal to the Hubble radius, 
Rh, centered at point x. We define 



/(A, r ) = inf \R 



J (Tw(\) — Tq°^ Xr,t ,x 



< 



rp(0) 

I 00 XR,t ,x 



,Vx eH 



111) 



In other words / is the smallest radius such that averaging over a ball of radius / centered at 
any point x lying within the Hubble volume relative to x always yields \5p\/p^ < 1. Note 
that our definition of Z(A,r ) depends on our choice of "window function" Xr,x ( x )- 

We now claim that /(A, r ) — > as A — > (and thus, in particular, / is always finite at 
sufficiently small A). To prove this, we note that if this result did not hold, we could find 
an Iq > and a sequence {A„ : n G N} converging to zero, such that l(X n , tq) > Iq for all 
n G N. Consequently, there would exist a sequence of points {x n } C T-L such that 



j ^00 (An) — Tqq^J Xlo,TQ,X n > j T^Xlo,TO,Xn 



(112) 



Since T-L is compact, there exists a subsequence — which we also denote as {x n } — converging 
to some z G %. By the triangle inequality, we have 



J (Too(K) - xi ,to,z + J (Tqo(K) - T$ j [Xh hn> ,x n 



Xlo,TQ,2 



> 



I T i0 K 

I ^00 XIo,to,x„ ■ 

(113) 

Taking the limit as n — > oo, we see that the first term on the left side vanishes because 
Too (A) converges weakly to Tqq and the second term vanishes by the lemma of section II. 
On the other hand, the right side is bounded away from 0, thus yielding a contradiction, 
thereby proving the desired result that /(A, r ) — > as A — > 0. 

Returning to the question posed four paragraphs above, since Too (A) provides the dom- 
inant contribution to the source term S ab , it seems clear that if, at finite A, V(A) is chosen 
to be so small that it does not include all source contributions lying within a homogeneity 
scale /(A, t ) about point x, then we cannot expect the source contributions from outside of 
V(A) to consistently average to zero to a good approximation. On the other hand, if V(A) 
is of order of the homogeneity scale or larger, then it seems plausible that (with suitable 
additional assumptions) a good approximation to t)( s ^ a P(\,x) will be obtained. Thus, we 
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have argued that fy( s ' a P(X,x) should be well approximated by 

l)W*(\,x) =4 J G^^M (V V (A,*') -Aw ; M^Mj fff^MdV, 

V(A) (H4) 
where V(A) may be taken to be a "cylinder" centered at x of proper spatial radius L and 
proper time "height" 2L, where L > 1(X,Tq). Note that in order to obtain the conclusion 
that we need only integrate over the local region V(A), it is essential that we have removed 
the "long wavelength part" of f) Q/3 . 

We now assume that /(A, To) <C Rc so that we can choose L such that L <C Rc, where Rc 
denotes the length scale of curvature of the background metric gfj (i.e., the Hubble radius, 
assuming spatial curvature is negligible). In the actual universe at the present time, we have 
Rc ~ 3000 Mpc and l(X, r ) < 100 Mpc, so we should easily satisfy all required criteria with 
L ~ 100 Mpc. In that case, it should suffice to take only the leading order terms in the 
Hadamard expansion for U , v , and V , v , as well as the leading order approximation to a 
in the Green's function expression ( lliop . This yields 



^ret ii'v 1 Vl U > 1 1 ^ I ~ V V 



a(r)r' 

+V^ lu/ (r,0,T,0)e(-T , + T-a(ry) , (115) 

where we have now put the field evaluation point x at the spatial origin of our coordinate 
system fll02p . Since V a ^ , v ,(t, 0, r, 0) is proportional to the curvature of g^, it is clear that 
the contribution of the second term will be down by a factor of (L/Rc) 2 from the first term, 
so we neglect this contribution. We also neglect "retardation effects", i.e., the difference 
between evaluating the source at time r — a(r)r' and time r. Our formula (I114p for (j^)"/ 3 
then reduces to 

h^(A,r,0)^4 J d^y^ 7 1 (s^(X, r, x') - A w 7 lim S ^^, x') \ 

(116) 

Next, we assume that Z(A,r ) not only goes to zero as A — > (as we have proven above) 
but is O(A) as A — > 0, i.e., /(A, t )/A remains bounded^ as A — > 0. Thus, if we choose L such 
that L/X remains bounded as A — > 0, by inspection of (11161) . it is then clear that terms in S ab 
that are o(l/A) as A — > will make only o(A) contributions to f)( s > a P. However, we assumed 



9 This precludes behavior wherein, e.g., Xoo(A) is 0(1) as A — > but its scale of spatial variation goes as 
A 1 / 2 rather than A. Such behavior would not be excluded by our previous assumptions. 
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above that as A 0, we have T 00 (A) = 0(1/A), T 0i (A) = O^/A 1 / 2 ) an d ^(X) = 0(1). We 
now make a final, additional assumption that terms of the form V c Vdf) afe are 0(1/A). In 
that case, it follows immediately that t ab is 0(1) as A — > (see ( 199]) ). It then follows that 
to obtain \)^ a $ to 0(A) accuracy, we need only take into account the contribution to S 
from Too- Thus, to 0(A) accuracy, we obtain 

[)^ 00 (A, r, 0) = 4 J dtf J ~ ] ^ (T 00 (A, r, x') - T^ 00 (r, x') - XT^ 00 (r, x')) a 2 (r)r' 2 dr' 

(117) 

and fjt 5 )^ = for all other components of f)( 5 ) a6 . We write 

<P = ~/ S)0 °. (118) 

Then <p{x) differs from the familiar formula for the gravitational potential arising in ordinary 
Newtonian gravity due to the matter lying within a ball of proper radius L about x only in 
the following ways: (a) There is a factor of a 2 (r) in (11171) . which arises from the trivial scaling 
difference between the spatial coordinates of (I102p and ordinary Cartesian coordinates, (b) 
2^(0)00 j g subtracted in the integrand of (11171) because ;he FLRW time slicing differs from 
a locally Minkowskian time slicing (see section IA of [l|); equivalently, the effects of T^ 00 
have already been taken into account via the dynamics of the FLRW background, (c) 
^j>(i)oo j g subtracted in the integrand of (11171) because its effects were already taken into 
account by 7^ ■ As discussed above, this subtraction of AT^ 00 gives the integral much 
better convergence properties. Thus, we conclude that the leading order short wavelength 
deviation from the FLRW background g^jj is described by Newtonian gravity, taking into 
account only the matter distribution lying within a region about x whose size is of order the 
homogeneity lengthscale. 

The motion of matter is given by 

= V a (\)T ab (\) = V a T a \\) + C a ac {\)T c \\) + C b ac (\)T ac (\) . (119) 

We may write 

C\ C (\) = C^\ C (X) + \C^\ C , (120) 
where, to leading order in wave map gauge, we have 

O (5)0 0l = V^, (121) 

C {s \o = VV, (122) 

O^^ = ^ ) V>-25 i o .V, ) 0, (123) 
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(with other components zero), and 

= l g[0)cd ^(J) + Vfe7 ( ? _ ^^L) j (124) 

The dominant terms in fll 19j) arise from C^- ) * 00 (A)T 00 (A) and correspond to the ordinary 
Newtonian gravitational effects on the motion of matter. Although, for small A, the contri- 
butions from C^ a bc will be much smaller than those arising from 0, it is important not to 
discard the terms in C^ a bc since they can produce large scale, coherent motions. 

C. Behavior of 7^ with dust source 

In this subsection, we will simplify the rather complicated equations for 7^ derived in 
section III under the assumption that l)^ ab (X) is of "Newtonian form". More precisely, we 
assume that the following quantities are uniformly bounded as A — > 0: 

{Vat) iS)0j , >&Vit) iS)0j , (125) 

The remaining components (f)^ *, f)*- 5 -**- 7 , and Vof) J ') are assumed to be o(A). These 
assumptions can be justified by generalizations of the arguments made in the previous sub- 
section, but here we will simply assume that they are valid. 

We will also assume that the matter stress-energy takes the form of "dust" 

T ab (X) = p(AK(AK(A) , (126) 

where u a (\) has norm —1 with respect to the metric 5w>(A), and p(A) > 0. (Recall that we 
have incorporated a cosmological constant into Einstein's equation, so the possible presence 
of "dark energy" of that form has already been taken into account. This assumption of the 
dust form of the stress-energy tensor as opposed to a more general form of non-relativistic 
matter is made here mainly for the purpose of obtaining definite equations involving familiar 
quantities.) We further assume that as A — > 0, u a (\) converges uniformly to u^ a , where in 
the coordinates of (11021) . we have = (1, 0, 0, 0). Note that since u a (X) — > uniformly 
as A — )■ 0, it follows by our lemma of section II and the positivity of p(A) that 

Ti? = P (0 M°M 0) , (127) 
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where p(°> = w-lim A _^ p(A). 

Let v a (X) denote the peculiar velocity of the dust relative to the "Hubble flow" u^ a , i.e., 
v a (X) is the projection of u a (X) orthogonal to u^ a in the metric g^. In accord with usual 
Newtonian limits, we assume that t> a (A)/A 1//2 is uniformly bounded as A — > 0. It follows that 

u a (X) =(l + h cd (X> { ° )c u {0)d + ^v c (X)v c (X)^j u^ a + v a (X) + o(A) , (128) 

or, equivalently, 

u a (X) = g ab (X)u b (X) = (l + h cd (X)u^ c u^ d + ±v c (X)v c (X)^j u ^+v a (X)+h ab (X)u^ b +o^) ■ 

(129) 

The results we shall now derive in this subsection will be rigorous consequences of the 
assumptions that have been stated above. We shall first show that, under the above as- 
sumption^, the quantities fi a bcdef, ^abclefghi an d v abcde a ^ van ish, so all of the "backreaction 

tensors" that appear in the equations for 7^ of section III vanish, except for /J^bcdef anc ^ 
K (i) 

^abcd- 

Taking the weak limit of \f d times (1981) (or, equivalently, using ( 1281) directly) we find that 
our Newtonian assumptions imply 



w-lim^09> = O, (130) 

A— >0 

where <p was defined by ( 11181) . Since the spatial metric is positive definite, this im- 



(S) 
ab 



plies that for any test function /, we have ||<9j(/0) ||£2 — > as A — > 0. Now, h 
^a)? — 5^06 f)^ C c + so ^e only possible contributions to u^cde fgh come from terms 

proportional to 

w-lim — 4>di4>di(j) . (131) 

A^O A 

Let / be any test function, and let g be a non-negative test function which is equal to 1 on 



10 In fact, we can show, without any Newtonian assumptions or gauge choice, that if tf^ = 0, then a a bcdef = 

and (^abodefgh = 0< However, this is not enough to tell us anything about v^bcde or the other components 
r (i) 

Ot Aabcdef Or UJ abcdefgh . 
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the support of /. Then we have 



w-lim — <; 

A^O A 



/d 4 x 



lim / —d>di(j)dj(j)f d x 
a^o / A 



lim 

A^O / A 



<Clim / |^(#)^(#)|d 4 x 
< Clim \\di(g<f>) \\ L 2 11^(^)11^2 

A — ^0 

= o. 



(132) 



Thus, we have ^bcdefgh = 0- Similar arguments show that fi a bcdef = 0. Finally, v^ilcde can 
only depend on terms proportional to 



w-lim — < 

A^O A 



w-lim di . 

a^o \ 2A 



(133) 



so it vanishes as well, as we desired to show. With this simplification, ( 1871) reduces to 



b)d 



lab 



87rT (1) 4- rv^ c d -?ttk^ c - SttkW c 4- ?7rn ( - 0) ( d - K^ cd \ 

° 7ll ab + " abed Z7FK ab c 87rK (a b)c + Z7[ 9 ab 1> cd K cdf 



(134) 



where 7^ = 7^ — hg^L l^ c C i an d where we have used the wave map gauge condition, 



V t£? = (see (HMD). 



Next, we evaluate K^ bcd - With our above assumptions of slowly moving dust matter, we 



have 



«2L = wJm i^ ) (A)p(A) Ma (A)^(A) = w-lim ^i?(A)p(A) . (135) 



The second equality follows from our lemma of section II because u c (X) — h^\X) /X — > 
uniformly as A — )■ and p(A) > 0. Using the Newtonian assumption that all components, 
except for the time-time component, of ^f b \X)/X converge uniformly to zero as A — > 0, the 
lemma of section II tells us furthermore that the only non- vanishing components of K^Lf are 



'-oooo 

-(1) 



-2 w-lim T p(X)(j)(X) 

A^0 A 

,(°) 



^oo = -2^w-hm-p(A)0(A) 



(136) 
(137) 
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To evaluate jJ^bcdef^ it * s more convenient to work with the quantity 



HVJ** = w-lim { V (a f^V 6) f^ - w 7 lim [V {a ^ cd V b) ^] 

A->0 A L A'^0 

= w-lim 4 \V a t) {s)cd V b t) {s)ef ] , 



(138) 



where the second line follows by the same type of argument as used above to show the 



ab 



29ab V c + O(E)f)), it is not difficult to see that 



vanishing of uj { ^ cdefgh . Since h$ 
ftabcdef can De expressed straightforwardly in terms of ^bcdef an< ^ vice- versa. Our Newtonian 
assumptions (1125)1 directly imply that the only potentially non-zero components of p^bcdef 



arc 



-(i) -(i) -(i) -(i) 

/^OOOOOO > /^OiOOOO ' /%000j ' ^ijOOOO 

-(1) -(1) -(1) 



(139) 



together with the components related to these by symmetries. However, it also follows that 
/^iiOOfcZ = °> since 



ffin/ /-V^gV.^d^ 



< Clim 

A^O 



J 00 



d 4 x 



< C'limllVK/^H^ 

A — 

= 0. 



(140) 



where in the second line we used that fact that Vjljjy /A is uniformly bounded as A — >■ 0. 

To further simplify p^b c defi we now a PP ea l to ( 1861) . With the simplifications arising from 
the vanishing of /J, a bcdef, u abode fgh' anc ^ ^afcLe together with the wave map gauge condition 
VJ) {s)ab = (see (ESD), we obtain 



±.-.(l)e -A ffK W _27rr; (0) KW £ 



From the form we have derived above for «Q it follows immediately that 



;i4i) 



S (l)< 



,-,(1)0 

f- aOjOfc 



0. 



(142) 



Again by the positive definiteness of the spatial metric, it follows that || Vj(/f) ^)/A 1 / 2 || L 2 — > 
as A — » 0, which implies /^jofco; = as well. By using similar Schwartz inequal 
arguments, it also follows that /^jooofc = ®- The wave map gauge condition V a f)^^ 



)o6 
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then yields 



POiOOOj 


~~ P- iOfeOj 


= 0, 


A^OiOOOO 


~~ h 1 iOjOO 


= o, 


Moooooo 


- u Wi 

— P OOiOO 


= 0. 



(143) 

The only nonvanishing components which remain are 

/4joooo = , (144) 

where 

En = w-lim \Vi(bVj(t) . (145) 

J A^O A 

Furthermore, by (1141ft and our previous expression for K^M, we obtain 

S* = -4vr w-lim -p(A)0(A) . (146) 

A^O A 

Finally, expressing ^i&Lfe/ m terms of p^bcdeji we fi nc l that the nonvanishing components of 

(i) 

Vabcdef are 

/4j0000 = ^^ij , (147) 

/Coo = 4^?^ , (148) 
„(l) _ 4 (0) (0)-.. (U q\ 

Next, we consider T a ^\ We have 



1 2 = tSpU t -m- t > 



(0) 
06 



wdimi(p(A)n a (A)n b (A)-p(°)ni )^ 0) 

A^0 A V 

w-lim ~ ( [p(X) - p (0) ] + p(X)v a (X)v b (X) + 2p(A) M [%(A) 

+ P (\)u^uf ) [MA)« (0)c « (p)d + « e (A)v c (A)] + 2 P (A) M (°)S[° ) V(A) 

+ Pff + 2«jM) ) + ^uf^ d cd + P^\f\^u^ 

p^uTnf + 2u$ P M +p W + u (o) u (o) p d)c c + u^uf^ + 2««V)c 

+M (0) M (0) T( )cd -( ? + 2T (0)c (a ^) _ 1 T (J). (150) 
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Here, we have introduced the quantities 

p« = w-lim [p(A) - p(0)} /X 

A->0 

= w-lWAK(A)/A, 



A->0 



and 



pS = w-Hm/)(A)f a (A)f6(A)/A 



A^O 



(151) 
(152) 

(153) 



Now, for a pressureless fluid, pu a is a conserved current, whose integrated flux through 
a spacelike hypersurface can be interpreted as (proportional to) the "number of particles" 
(often referred to as the "number of baryons") in the fluid. The perturbation, of the 
density (relative to the background metric gjj£) of particles on a r = const, hypersurface is 
given by 



_(i)Jo) _ 1 
PM^ c = w-hm- 



T de (X)u d (XX(X)e abc (X) - TjJV°> V> e e2>" 



(154) 



where n a (X) and e abc {X) are the unit normal and 3- volume of the hypersurface in the metric 
g a b{X) and n^ a and e ® b \ are the corresponding quantities for g^ b \ We obtain 



(155) 



We now return to our perturbation equation (11341) for 7^. We use the above explicit 
expressions for fi^lL anc ^ /^i&Lfe/' we use (H5Up to substitute for and we use (j!55p to 
eliminate p^ in favor of p^r • Bringing all terms explicitly involving to the left side and 
using the background Einstein equation, we obtain 



(a7fo) c 



! a6 7 V 



EE 87T0, 



where the components of Q ab are explicitly given by 

ft - n (1) 4- 1 ti( 1 ) < 1 



e 



Oi 



P; 



(1) 



1 - _ J_JO)^k 



(156) 

(157) 
(158) 
(159) 



47T ' 87T 

This expression for the "effective perturbed stress-energy" of matter agrees with the expres- 



sion obtained by 



28| (see also |2l|) 
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All of the terms on the left side of (11561) together with the terms p M and P!> appearing 
in <d a b would be present in ordinary linearized perturbation theory. We previously showed in 
section II that, in the absence of gravitational radiation, small scale inhomogenities cannot 
affect the dynamics of the background metric. We now see that, under the assumptions 
stated at the beginning of this subsection, the only effect that small scale inhomogeneities 
have on long wavelength perturbations is to add the terms involving and Hjj to the 
effective perturbed stress-energy Q ab . These additional terms provide precisely the contri- 
butions to the energy density and stresses that one would expect from kinetic motions and 
Newtonian gravitational potential energy and stresses. In particular, for the energy den- 
sity, these terms have the effect of shifting the proper mass density appearing in the FLRW 
background to an "ADM mass density" (to first order in A). Although the presence of these 
terms is very important as a matter of principle, they should be extremely small compared 
with p$ (and even with ). 

Finally, we note that although we have written ( 11 56ft in wave map gauge (I104p . we may 
use the gauge freedom with respect to smooth, one-parameter groups of diffeomorphisms 
discussed in the paragraph below (157)1 to impose any other desired gauge condition on 7^ . 
Since the left side of ( I156P together with the terms p$ and in Q a b have the same gauge 
transformation properties as in ordinary linearized perturbation theory, these terms will take 
the same form as in ordinary linearized perturbation theory when transformed to the new 
gauge. On the other hand, the terms in Q ab involving and are gauge invariant under 
the transformations induced by a smooth, one-parameter group of diffeomorphisms. This 
enables one to write down the form of fl 1 5 6 [) in other gauges. 



V. APPLICABILITY OF OUR FORMALISM TO THE REAL UNIVERSE 

In this paper, we have developed a formalism that enables one to take full account of small 
scale non-linearities in Einstein's equation. This formalism is, in essence, an adaptation of 
Burnett's formulation of the "shortwave approximation" for analyzing the self-gravitating 
effects of short- wavelength gravitational radiation. The key idea of our formalism is to 
consider an idealized, one-parameter family of metrics, g a b(A), with the property that, as 
A — > 0, the deviation of g a b(ty from a "background metric" goes to zero in proportion to 
A, but the scales over which the metric and stress-energy vary also go to zero in proportion 
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to A. In principle, this allows the small-scale nonlinear terms in Einstein's equation to have 
a significant effect on the large scale dynamics of gf b \ However, we proved in section II that 
the only such effects that can actually arise correspond to the presence of an effective stress- 
energy associated with gravitational radiation. No new effects can arise from the presence 
of matter (provided that the matter satisfies the weak energy condition), and no effects can 
arise that in any way mimic "dark energy" . 

In sections III and IV, we applied our formalism to analyze the leading order deviations 
of g a b(X) from g^} . Within our formalism, these deviations can naturally be split into 
"long wavelength" and "short wavelength" parts. The long wavelength part satisfies a 
modified version of the usual linearized perturbation equation. We argued that, in the case 
of cosmological perturbations with nearly Newtonian sources, the short wavelength part of 
the leading order deviation of g a b{X) from g^ is described by Newtonian gravity, taking into 
account only the local matter distribution. 

However, our actual universe is not an idealized limit of spacetimes with inhomogeneities 
on arbitrarily small scales, but is a particular spacetime with finite amplitude, finite wave- 
length deviations from a FLRW model. Quantities in our formalism like the "long wavelength 
perturbation," 7^, are defined by taking weak limits as A — > 0. How do quantities like 7^ 
that arise in our formalism correspond to quantities observed in the actual universe? Fur- 
thermore, to what extent are we justified in applying our results to the actual universe, i.e., 
how do we know whether the actual universe is sufficiently "close" to the idealized limit 
considered in our formalism that results derived using our formalism should hold to a good 
approximation? 

The above questions can also be asked in the context of ordinary linearized perturbation 
theory. In this context, one also deals with a limit of a one-parameter family of metrics 
g a b{X, x) (now assumed to be jointly smooth in (A, x)) and the formalism obtains results for 
idealized quantities like dg a b(X, x)/d\\\ =0 . Nevertheless, in the case of ordinary perturbation 
theory, the above questions can be answered in a relatively straightforward manner. In an 
actual ("finite A") spacetime with metric g a b, one may introduce a background metric g^J and 
identify the difference, h a b = g a b~9^b '> w hh Xdg a b(X,x)/dX\\ =0 . A fundamental criterion for 
the validity of ordinary perturbation theory is that <C 1 in some orthonormal basis of 
gf b \ However, we also need to satisfy conditions that state that first and second spacetime 
derivatives of h^ v are sufficiently "small," since these quantities appear in the nonlinear 
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terms that are being neglected in Einstein's equation. In particular, the criteria for the 
applicability of ordinary linear perturbation theory cannot be satisfied in a cosmological 
spacetime with \5p\ > po, as occurs in the real universe. Indeed, it is for this reason that 
we have developed the formalism of this paper. It should be noted that even in the case of 
spacetimes that do satisfy the basic criteria for the validity of ordinary linear perturbation 
theory, it would be very difficult to obtain precise error estimates for this approximation. 

If one wishes to apply the formalism of this paper to a particular spacetime, such as 
our universe, it also is necessary to introduce a background metric . As in ordinary 
perturbation theory, it is essential that the difference, h ab = g ab — gf b \ satisfy \h^ v \ <C 1 
in some orthonormal basis of g® b . The main advantage of our formalism over ordinary 
perturbation theory is that it imposes much weaker restrictions on spacetime derivatives 
of h^ u . In particular, for cosmological spacetimes, it allows \8p\ ^> p 011 small scales. 
However, it is clear that some further restrictions in addition to j/t^l <C 1 must hold for 
the results derived using our formalism to be a good approximation. Although it would be 
extremely difficult to formulate mathematically precise criteria for the validity of applying 
our formalism to a given spacetime — much more difficult than for ordinary perturbation 
theory — we now shall propose a rough criterion for its validity. 

In essence, taking the weak limit as A — > of quantities that arise in our formalism 
corresponds to taking spacetime averages of these quantities over arbitrarily small regions 
of spacetime before letting A — > 0. For a given, fixed spacetime, such as our universe, the 
weak limits in our formalism should thus be identified with spacetime averages over small 
regions. Thus, for example, the quantity Xj^ (x) [where 7^ = w-lim^o ^o&(A)/A] should 
be identified with the spacetime average of h ab over a suitable region, 1Z X , centered on x. 
In order that the spacetime under consideration be "sufficiently close" to the idealized limit 
A —7- of our formalism that our results should apply, it is necessary that 1Z X be small 
compared with the curvature scale, Re, of gj^; indeed, the notion of "averaging" would be 
highly ambiguous if this were not the case. However, it also is necessary that TZ X be large 
enough that our perturbative approximations should apply. In particular, X'y^ must satisfy 
the same "smallness" conditions as required in ordinary linearized perturbation theory, i.e., 
its spacetime derivatives must be appropriately small. Similarly, the perturbative quantity 
ATqq- 1 (see (I7H1) ) should correspond to the spacetime average of T 00 — TqqK In order that our 
approximations should apply, we must choose 1Z X large enough that |ATq^| <C . 
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The condition that the spacetime average of T 00 — T q be less than T q was considered 
in section IVB above and used to define the homogeneity scale, / (see fillip V Similar 
"homogeneity scales" can be defined for all other "long wavelength" quantities arising in 
the analysis of section III. Since our averaging region must be large compared with these 
homogeneity scales but small compared with the scale of the curvature of the background 
metric, it is clear that a necessary condition for the applicability of our formalism is that 
these homogeneity scales be small compared with the curvature of the background metric. 
We believe that this condition, together with <C 1, should also be sufficient for the 

applicability of our formalism, since it should guarantee the appropriate "smallness" of all 
long wavelength perturbative quantities. 

In the case of our universe, it seems reasonable to assume that the only relevant homo- 
geneity scale is the density homogeneity scale, I, defined by fillip . In the present universe, 
we have / < 100 Mpc whereas the scale of the curvature of the background metric (i.e., the 
present Hubble radius) is of order Rq ~ 3000 Mpc. Thus I <C Rc- In addition, common 
sense estimates indicate that, except in the immediate vicinity of black holes or neutron 
stars, we have {h^l < 10 -5 . Thus, if we define "long wavelength perturbations" by averag- 
ing on scales of order L ~ 100 Mpc, the criteria we have proposed for the validity of our 
approximations should be satisfied. Consequently, we believe that the results of this paper 
should be applicable to our universe to an excellent approximation. 

If this is the case, then the following conclusions can be drawn: (1) As shown in section 
II, the only effect that small scale inhomogeneities can have on the leading (i.e., zeroth) 
order large scale dynamics of our universe is that of a P = |p fluid, corresponding to the 
presence of gravitational radiation. In particular, small scale inhomogeneities cannot mimic 
the effects of "dark energy." (2) The deviation, h a b, of the metric from a FLRW model can 
be broken up into "long wavelength" and "short wavelength" parts. The long wavelength 
part corresponds in our universe to averaging h a b over a spatial scale of order L ~ 100 Mpc. 
As analyzed fully in section III and subsection IVC, it satisfies a linear equation with an 
additional source term due to the short wavelength part. As argued in subsection IVB, the 
short wavelength part should be described by Newtonian gravity, taking into account only 
the matter within a proper distance of order L ~ 100 Mpc of the point under consideration. 

Thus, the analysis of this paper goes a long way towards justifying many of the key as- 
sumptions made in cosmology. In particular, if the matter in the universe is non-relativistic, 
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it suggests that to calculate structure formation in the universe on all scales — from the Hub- 
ble radius or larger down to arbitrarily small scales — it should be a good approximation to 
evolve the long wavelength part of the deviation of the metric from a FLRW model by ( 11561) . 
to calculate the short wavelength part by ( I117|) . and to calculate the motion of the matter 
by (HUD together with (THUjl -dTm 

In addition, our analysis suggests how further improvements can be made to get more 
accurate approximations. To improve upon the description of the "long wavelength" part of 
the deviation of the metric from a FLRW model one could go to higher order perturbation 
theory. In particular, the second order correction to 7^ would be defined in our framework 
by 

/u (0) x (L) 
(2) = w _ lim flab (A) - 9ab ~A 7afe (16Q) 

a0 A^O A 2 

(2) 

assuming, of course, that this weak limit exists. The equations satisfied by 7^ — taking 
full account of the small scale inhomogeneities — could be derived by the methods used in 
section III. The description of the short wavelength part of the deviation of the metric from 
a FLRW model could be improved by a more accurate treatment of the integral relation 
( 11081) . However, the investigation of such improvements is beyond the scope of this paper. 
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